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EDITORIAL 


Como se pode constatar pela simples observação das capas das últi- 
mas revistas esta estava a ser publicada com dois meses de atraso. 


Tentamos divulgar as razões que levaram a este estado de coisas no 
último Editorial e como se viu residem sobretudo nas dificuldades havidas 
no arranque das nossas actividades após longo tempo de encerramento 
compulsivo. 


Desde artigos prontos a publicar e que ficaram ultrapassados a originais 
perdidos, a anunciantes que foi necessário contactar de novo até ao aumento 
substancial das matérias primas e mão de obra agravado ainda pela não 
totalmente solucionada crise financeira da AEIST tudo foram problemas 
que tentamos e conseguimos em parte superar. 


Neste momento com a maioria dos problemas ultrapassados e com a 
secção em funcionamento tentávamos a toco o custo recuperar os meses 
em atraso e pensávamos em Julho próximo publicar a revista referente a 
esse mesmo mês. 


Verificámos, no entanto, que foi muito grande o número de revistas 
devolvidas com multa e tivemos conhecimento que muitos assinantes tiveram 
a amabilidade de a pagar, aos quais pedimos desculpa por isso. 


Contactados os CTT fomos informados que só poderiamos beneficiar 
das facilidades por eles concedidas no envio de publicações periódicas 
sempre que elas sejam enviadas com um atraso máximo de um mês. Esse 
atraso é medido em relação ao mês a que a publicação se refere conforme 
indicação na capa. 


Não nos sendo possivel ultrapassar mais esta dificuldade a acrescentar 
ao tremendo esforço que se estava a dispender para tentar recuperar deci- 
dimos saltar os meses de Dezembro e Janeiro sendo esta revista pois refe- 
rente a Fevereiro(/5). 


Ao tomar esta atitude tivemos presente que em nada vamos prejudi- 
car os Assinantes na medida em que na assinatura é especificado o número 
de exemplares a que ela se refere e não qualquer período de tempo. 


Agradecemos uma vez mais a vossa compreensão. 


A Técnica 
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BASES DE ESTUDO ECONÓMICO DE UM 
EMPREENDIMENTO MINEIRO 


RESUMO 


A análise económica aplicada a empreendimentos 
mineiros pressupõe o domínio de uma sequência de 
conhecimentos interligados, respeitante a informação, 
avaliação e decisão. A validade daquela análise é lógi- 
camente condicionada pela representabilidade destas 
actividades. Neste trabalho, faz-se referência à neces- 
sidade de obtenção de dados de confiança, mensionam- 
-Se os vários métodos deterministicos e probabilísticos de 
avaliação económica e apontam-se critérios de decisão 
mais utilizados, com especial realce para os que se 
apoiam na simulação por computador. 

A natureza interdisciplinar deste estudo é posta em 
relevo, embora certos aspectos da análise sejam típicos 
da Indústria Mineira. 


| — INTRODUÇÃO 


O estudo económico de um empreendimento mi- 
neiro é uma análise que depende de vários factores, uns 
de natureza técnica, outros de natureza financeira, 
outros até de natureza comercial (mercados, por 
exemplo). 

Um estudo económico bem conduzido deve incluir 
o custo inicial no empreendimento, assim como os lucros 
esperados no futuro, os quais são obtidos a partir de 
uma análise do fluxo de capitais entrados e saídos na 
Empresa durante a vida prevista para o empreendimento. 
Esse fluxo entrará em consideração, por um lado, com 
as receitas brutas da exploração, por outro, com os 
custos operacionais e os custos indirectos (tais como 
depreciação, amortizações, deplecção do jazigo, custos 
de aquisição do capital, etc.). 

A análiss econômica do empreendimento faz uso 
dessas grandezas, calculadas a priori, e inclui as mes- 
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C. DINIS DA GAMA 
Prof. Catedrático 
da Universidade de Luanda 


SUMMARY 


The economic analysis of mining prospects is 
conducted through the knowledge of a series of topics 
concerning information, evaluation and decision. The 
validity of each analysis depends upon the reliability of 
those forecasting activities. 

In this paper, the need for obtaining reliable data, 
as well as for applying the available methods, both 
deterministic and probabilistic, of economic appraisal, 
and also to utilize decision criteria are underlined. Spe- 
cial emphasis is ascribed to computer simulation 
methods. 

The interdisciplinary nature of this economic ana- 
lysis is mentioned, although some aspects are chara- 
cterístic of the Mining Industry. 


mas em métodos econômicos já consagrados para o 
estudo da rentabilidade de investimentos. Através do 
emprego destes métodos, encaram-se os diferentes graus 
de interesse que caracterizam as várias alternativas dis- 
poníveis para conduzir um mesmo empreendimento. 

O risco inerente a cada uma das alternativas tam- 
bém pode ser quantificado, permitindo assim a escolha 
da solução que minimize a probalidade de insucesso 
tou risco), simultaneamente com o lucro máximo que 
se obterá potencialmente através das várias alterna- 
tivas, 

A melhor alternativa será normalmente aquela que, 
entre todas as possiveis, maximizar o lucro previsto para 
o empreendimento, dentro de um minimo de risco. 

À incerteza dos resultados é consequência da limi- 
tada informação que habitualmente existe na indústria 
mineira, no que diz respeito ao cálculo das reservas, ao 
conhecimento do jazigo a explorar, aos custos de explo- 
ração previstos e à evolução futura dos mercados. O 
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carácter aleatório destes números implica que se utili- 
zem métodos de análise probabilística apropriados para 
o respectivo processamento. 

Ântes de passar em revista esses métodos, come- 
car-se-ãá por uma breve descrição dos métodos conven- 
cionais de estudo da rentabilidade de empreendimentos, 


2 — ESTUDO DETERMINÍSTICO DA EXPLORABILIDADE 


Um empreendimento mineiro inicia-se geralmente 
com a realização de um conjunto de trabalhos de pros- 
pecção a pesquisa, ou com a aquisição de dados geoló- 
gicos obtidos por outrém, respeitantes a determinada 
área mineralizada, 

No termo de essa fase preliminar de prospecção 
(maior ou menor, conforme o capital disponivel) põe-se 
invariavelmente um problema de decisão sobra a conti- 
nuação ou paralização dos trabalhos. Como para rea- 
lizar o reconhecimento e valorização do jazigo são indis- 
pensáveis maiores investimentos de capital, é necessá- 
rio assegurar que os mesmos possam vir a ser recupe- 
rados através de uma futura exploração. 

Este estudo envolve, portanto, uma análise sobre a 
viabilidade económica do empreendimento que, se for 
favorável, prossegue com a selecção da melhor das 
alternativas que se deparam à continuação das ope- 
rações. Mais ainda, urge demonstrar que a alternativa 
escolhida é mais atractiva para os investidores que 
outras actividades geradoras de dinheiro com menores 
riscos. 

Frequentemente, os cálculos mineiros são elabo- 
rados em função de dois parâmetros que se podem 
prever através dos trabalhos de prospecção e pesquisa. 
São eles: 


— a produção prevista para a futura mina; e 
— q teor recuperável do minério, 


Com base na informação existente e na experiên- 
cia, fazem-se estudos de rentabilidade, avaliando as 
futuras receitas em confronto com as estimativas de 
custos, empregando métodos típicos de Análise Eco- 
nómica. 

Vejamos quais os méodos mais correntes que se 
usam para a avaliação das oportunidades de investi- 
mentos: 


- Método do periodo de reembolso; 
— Método da taxa de reembolso; 
— Método do valor actual; 

— Métodos de desconto. 


2.1 — Método do periodo de reembolso 
É conhecido pelo método do payback ou payout e 


consiste simplesmente em calcular o tempo que um 
projecto leva a gerar receitas iguais ao capital investido 


O = 


normalmente conhecido por inflacção. 
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para o pôr em execução, É pois um método que visa a 
determinação do número de anos que o projecto leva 
até repôr o seu custo inicial. 

Teoricamente, dentro dos princípios deste método, 
um projecto seria preferivel em relação a outro se fosse 
pago em menos anos, mas este critério é incompleto, 
pelas seguintes razões: 


— Não toma em consideração o valor actual do di- 
nheiro, ou seja, a influência do factor tempo no valor 
do capital. 

Como se sabe, o valor actual V, de certa quantia 
exprime-se em função do seu valor futuro V, ao fim 
de n anos, por: 


? Ve 
/ +" 


sendo | o juro anual, expresso em forma decimal (*). 

Por exemplo, o valor actual da quantia (ou do 
lucro) de 1000 contos realizada em oito anos, utilizando 
um juro de 10% ao ano, é: 


1000 


—— > = 466,5 contos 
(1 + 0.1) 

— Qutro inconveniente é o método não se interes- 
sar pelos lucros que se conseguem após o periodo de 
reembolso considerado, 

—— Também não serve como critério de explorabili- 
dade, como se deduz do seguinte exemplo: suponhamos 
duas alternativas A e B que exigem igual custo inicial. 
Se a À tiver um periodo de recuperação de três anos e 
uma vida produtiva de cinco, e B um período de recupe- 
ração de cinco anos e uma vida produtiva de vinte, se- 
gundo este método A seria preferível, embora fornecesse 
lucros muito inferiores no seu conjunto. 

Também pode acontecer que A e B tenham iguais 
investimentos e iguais períodos de reembolso, mas 
que um deles recupere mais depressa que o outro uma 
percentagem importante (B0% por exemplo) do inves- 
timento. Logicamente, esse seria preferivel, embora 
segundo este método estivessem em iguais circunstân- 
cias. 

Por estes motivos, o método do periodo de reem- 
bolso é considerado pouco próprio para análise de 
investimentos complexos, como são os da Indústria 
Mineira. 


2.2 —- Método da taxa de reembolso («rate of 
return» ) 


Consiste essencialmente em determinar o quociante 
entre o lucro anual médio e o capital inicial despendido. 
Esse quociente constitui a taxa de reembolso, que quanto 
mais alta for mais conveniente é a alternativa, 


[*) Este conceito não toma em consideração a diminuição com o tempo do poder aquisitivo do dinheiro, fenómeno inevitável, 


TÉCNICA 423 


PROJECTO 


Receitas 


A Custos 


Lucros 


Receitas 


B Custos 


Lucros 


TAXA DE 


VALOR ACTUAL 


Receitas 


C Custos 


Lucros 


Tal como o anterior, este método também tem o 
defeito de não distinguir entre projectos que obtenham o 
mesmo lucro total no mesmo intervalo de tempo, embora 
uns apresentam lucros crescentes outros decrescentes 
e outros ainda constantes, Para ilustrar esta afirmação, 
consideremos o exemplo de três projectos 4, Be £€ 
que apresentam a mesma vida útil (cinco anos), os 
mesmos lucros totais e o mesmo capital inicial (400 
contos, por exemplo). 

Enquanto o projecto 4 tem receitas constantes, B 
tem receitas sempre crescentes e € sempre decrescentes. 
Embora tenham todos a mesma taxa de reembolso, 
é fácil compreender que C é melhor que À pois recupe- 
ra mais cedo a mesma quantia. Pela mesma razão, A 
é preferivel a B. 

Embora estas conclusões sejam intuitivas, elas não 
fazem parte deste método, mas serão quantificadas pe- 
los métodos que se indicam a seguir. 

Outro aspecto que o método da taxa de reembolso 
não considera é o custo do dinheiro. De facto, se o 
juro a que foi adquirido é elevado (comparável aos 
259%), qualquer dos projectos não é rentável. Na reali- 
dade, entrando com esse critério noutros métodos, pode 
concluir-se que os projectos são rentáveis e se classi- 
ficam por ordem decrescente de interesse em CL, A e B, 

Concluímos pois que o método da taxa de reem- 
bolso satisfaz como critério de explorabilidade, por ser 
contra os projectos a curto prazo, ao contrário do mé- 
todo do período de reembolso, 


2.3 — Método do valor actual 


Já foi definido o conceito de valor actual, De 
facto, o dinheiro possui um valor que é função do 
tempo e a recuperação de investimentos ao fim de al- 
guns anos deve entrar em consideração com esse facto. 
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Assim, por exemplo, um lucro de 100 contos obti- 
do ao fim de um ano (após qualquer investimento) 
pode ser reinvestido e valer 150 contos ao fim de cinco 
anos. Por outras palavras, um lucro de 100 contos ao 
fim de um ano equivaleria ao lucro de 150 contos ao 
fim de cinco anos, supondo inexistente o problema 
da inflacção. 

O valor do dinheiro em diferentes intervalos de 
tempo depende, pois, de 


— Taxa de juro (normalmente escolhe-se um juro 
igual ao que onera a aquisição do capital); 


— Quantidade de dinheiro investido; 
— Tempo ao longo do qual é investido. 


Assim, usando a expressão V | = V,/(1 + i), esta- 
belece-se quanto vale em cada um dos anos seguintes 
a unidade monetária, e multiplica-se pelos lucros con- 
seguidos (ou previstos) ao longo desses anos. O soma- 
tório desses lucros, corrigidos para o seu valor actual, 
representa a rentabilidade do investimento, que pode 
antão ser comparada com o seu custo inicial. A dife- 
rença entre aquele somatório e o investimento (ambos 
referidos aos seus valores actuais) representa o valor 
actual líquido. 

Entre várias hipóteses de investimento, escolher- 
-se-ã, evidentemente, a que tiver maior valor actual 
líquido. 

Define-se também a taxa de valor actual como a 
relação entre o valor actual líquido e o investimento 
inicial, 

No caso de o investimento não ser intantâneo, mas 
se concretizar ao longo de p anos (anos de pré-produ- 
ção), a taxa de valor actual será o quociente 
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Valor actual Lucros, em valor actual, dos projectos 


ANOS de 1 c. à | 
taxa de 7% A | B c 

1 0,934 93,4 0,0 186,8 
2 0,873 87,3 43,0 130,9 
3 0,816 81,6 81,6 81,6 
q 0,762 76,2 114,3 38,1 
5 0,713 71,3 142,6 0,0 
SOMA 410,0 402,1 436,0 
INVESTIMENTO 400,0 400,0 400,0 

Valor actual líquido 10,0 2,1 36,4 | 

Taxa de valor actual 2,5% 0,5% 3,1 | 


| 


x=n V, benefício utilizável como dividendos, como capital para 
2 PRE outros investimentos, etc. 

x | Ao longo da vida de um empreendimento, é habi- 

=p Va tual ter no início um fluxo de capitais negativo (período 

Ph = de pré-produção) que passa a positivo após a produção 

x dedo, dar origem a receitas superiores aos custos reais (cus- 


em que n é o número de anos do empreendimento, | o 
juro anual pré-estabelecido e V, é o lucro futuro no 
ano x. 

Entre várias alternativas, será preferida a que 
apresentar maior taxas de valor actual, 

Retomemos o exemplo dos três investimentos A, 
B e €C, admitindo que o juro pré-fixado era de 7% ao 
ano. 

Observa-se, pois, que a alternativa mais conve- 
niente é C, seguida de A e finalmente de B. 


2.4 — Métodos de desconto 


Sob esta designação, incluímos os chamados meé- 
todos de fluxo de capitais descontados [discounted 
cash flow. 

O fluxo de capitais num empreendimento não é 
mais que o balanço contínuo entre as despesas, soma- 
das a outros custos exigidos pelo arranque desse em- 
preendimento, e as receitas de caixa que são esperadas. 
Tal balanço é, em geral, estabelecido ao ano. 

Nos custos de um empreendimento há que incluir 
as despesas de capital e de lançamento do projecto, 
enquanto que nas receitas se incluem lucros, deprecia- 
ções e a deplecção (*). 

Os métodos de desconto trabalham sobre a geração 
de dinheiro e não na contabilidade dos lucros, porque 
só quando o dinheiro é recebido é que constitui um 


tos operacionais + impostos + dispendios de capital). 

Ássim, em cada intervalo de tempo, um empreen- 
dimento pode ser caracterizado por fluxos de capital 
positivos ou negativos, conforme a soma algébrica das 
receitas e das despesas for positiva ou negativa, dentro 
desse mesmo periodo, 

Só se podem apreciar em conjunto esses saldos 
(positivos ou negativos), se os situarmos no mesmo refe- 
rencial de temipo, ou seja, se calcularmos os respectivos 
valores actuais. 

Designa-se por desconto a conversão de fluxos 
de capital futuros em valores actuais. 

Define-se taxa de desconto interna (ginternal rate of 
return») como o juro que aplicado aos saldos positivos 
lhes confere valores actuais correspondentes aos saldos 
negativos, quando descontados ao mesmo juro, 

Portanto, a taxa de desconto interna não é mais do 
que o juro que faz anular a soma algébrica de todos os 
valores actuais que caracterizam a vida de um empreen- 
dimento, 

Analticamente, a taxa de desconto interna r é dada 
pela relação 


n 
v Vx 


| tom 


— | 
x=p+1i (1+r) x 


onde as letras têm o mesmo significado já apontado a 
propósito do método do valor actual. 


[*) Considera-se deplecção a amortização do valor inicial do jazigo, atendendo ao carácter esgotávvel do mesmo. Exis- 


tem duas técnicas de cálculo de deplecção: por custos e por percentagens admissíveis do lucro bruto. 
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ES a 


Note-se que, se o investimento é instantâneo 
p= 0), o seu valor inicial é simplesmente dado por 


V V V, 
v a" Hi ra 12 n 


— + +... +——— 
T+) + (1+r)" 


Maior incerteza 
onda r traduz a taxa de desconto interna nesse caso. 
Como critério de selecção entre várias oportuni- 
dades de investimento, escolhe-se aquela que tiver 
maior taxa de desconto interna. 
Normalmente r é calculado por tentativas e in- 
terpolações, 


Menor incerteza 


PROBABILIDADE 


VALOR ESPERADO 


Vejamos o exemplo apresentado: Fig. 1 
CASO A CASO B CASO €C 
ANOS a pt 
r= 0% r=78 r= 10% r= 10% r=12% 
1 90,9 93,4 0,0 181,8 178,6 
2 82,6 87,3 43,0 123,9 119,5 
3 75,1 81,6 715,1 75,1 71,2 
4 68,3 76,2 102,0 34,1 31,7 
5 62,1 71,3 124,2 0,0 0,0 
379,0 410,0 344,3 414,9 401,0 
T. D; 1. | 
lada) 


Os diferentes projectos serão escalonados por 
ordem decrescente de valores da taxa de desconto in- 
terna, preferindo-se o caso C, depois 4 e finalmente B. 


3 — ESTUDO PROBABILÍSTICO DA EXPLORABILIDADE 


Os critérios da explorabilidade apresentados ante- 
riormenta são avaliados a partir da definição de pará- 
metros constantes (reservas, custos, receitas, lucros, 
ate.) que teoricamente correspondem a previsões rigo- 
rosas, 

Supondo que as grandezas desses parâmetros são 
valores médios (ou valores esperados), o critério ana- 
lítico mais não permite que achar um resultado espe- 
rado. 

No entanto, estas estimativas nunca podem ser 
rigorosas, pois em todas elas há uma incerteza maior ou 
menor. Na prática, as alternativas que se põem à explora- 
bilidade de um empreendimento mineiro medem-se não 
só em função dos lucros esperados, mas também do grau 
de incerteza que os caracteriza (ou risco). 

A incerteza é dificil de medir quantitativamente, 
sendo habitual defini-la a partir da distribuição de proba- 
bilidada em torno de um valor médio ou esperado. Usan- 
do a curva normal de Gauss, podem caracterizar-se parê- 
metros com maior ou menor grau de incerteza (Fig. 1). 
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As distribuições de probabilidade que se obtêm 
para caracterizar resultdos obtidos a partir do cálculo 
de vários parâmetros aleatórios, podem ser obtidas por 
dois métodos distintos: o método analítico e o método 
de Monte Carlo. 


3.1 — Método analítico 


Os parâmetros mais influentes dum estudo econó- 
mico podem ser caracterizados por uma estimativa do 
seu valor médio (ou valor esperado) e por uma distri- 
buição de probabilidade associada à respectiva variação. 

Admitindo a distribuição de probabilidade normal (o 
que está de acordo com numerosas determinações expe- 
rimentais), podem estabelecer-se intervalos de confiança 
das variações dos diferentes parâmetros em torno do 
seu valor médio, 

Tais intervalos são escolhidos arbitrariamente em 
função das previsões que se fazem, traduzindo-se por 
elevadas probabilidades (de BO, 90, 95%) para que os 
desvios em relação ao valor médio não ultrapassem 
uma certa função do desvio padrão, que é dada pelas 
áreas situadas debaixo da curva normal. 

Assim, por exemplo, existe uma probabilidade da 
68% para que os valores a, estejam compreendidos entre 
im — s) e tm + s) sendo m o valor médio e s o desvio 
padrão: 


115 


O intervalo de confiança dos 90% limita os valores 
antro (m — 1,645) e im + 1,645) enquanto que o in- 
tervalo de 95% estabelece os limites (m + 25), sendo 
estes dos mais utilizados. 

Vejamos um caso de aplicação ao estudo de um 
empreendimento mineiro. Se os custos esperados por 
tonelada forem os indicados na Tabela a seguir, apre- 
sentando-se ainda as suas variações dentro do intervalo 
de confiança de 90%, qual será o custo total e a res- 
pactiva margem de variação? 


| 


Custo 
OPERAÇÕES | MEpársido 
| S/ton 
Exploração mineira r 4500 
| Beneficiação | 1520 
| Transportes 0OSBO 
Custos fixos | 0575 
| Custos de administração | 0540 
| [—— 
| TOTAL | 6595 


O problema resolve-se simplesmente aplicando 
uma conhecida regra do cálculo de probabilidades (*), 
segundo a qual o desvio padrão do custo total é igual 
à raiz quadrada da soma dos quadrados dos desvios 
padrão dos custos componentes. 

será, portanto, 


/ 0861º + 0806” + 0503º + 0812º + 0806! = 0868 


Assim, o intervalo de confiança dos 90% para o 
custo total, é: 


Fá 5583 


6595 + 1,64 X 0$68 = 6$95 + 1812 d 
> 8807 


Existirá, pois, uma elevada confiança (definida por 
uma probabilidade da 90%) para que o custo total fique 
situado entre 55883 a 8507, 

Esto processo é suficiente jpara tratar problemas 
em que todas as variáveis têm distribuições normais, 
sendo, além disso, independentes umas das outras. 

Em casos mais complexos e onde é necessário 
obter a avaliação econômica de um processo completo 
(por exemplo, um empreendimento mineiro), o método 
analítico é difícil de aplicar. 


(*) O chamado teorema da adição das varâncias. 
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Recorre-se, na maioria dos casos complexas, ao mé- 
todo da Monte Carlo. 


3.2 — Simulação pelo método de Monte Carlo 


O método de Monte Carlo é um processo de amos- 
tragem onde se fazem aproximações realisticas de fun- 
ções aleatórias, tais como os critérios de explorabilidade. 

Desde que existam estudos bem conduzidos sobre 
os valores esperados e as variâncias dos parâmetros 
fundamentais do problema, pode fazer-se um conjunto 
muito grande de simulações do cálculo real, analisan- 
do-se posteriormente a distribuição de probabilidade dos 


Estimativa de variação 
dentro do intervalo de 
confiança de 90% 


Desvio padrão 
(3) = [2)/1,64 


+ 1500 0560 
= 0810 0506 
+ 0805 0503 
+ 0520 0512 
+ 0810 0506 
+— 1912 0568 


resultados, à custa do tratamento estatístico dessa con- 
junto de soluções. 

As distribuições de probabilidade de cada variável 
não precisam de ser normais, pois o que é importante 
é obter uma le! de probabilidade acumulada, de forma 
que a saída de qualquer número aleatório compreendido 
entre O e 1 corresponda a determinado valor dessa variá- 
vel. À custa da saída de números aleatórios para todas 
as variáveis, obtem-se o valor delas em cada caso e 
realizam-se os cálculos necessários para achar o critério 
da explorabilidade correspondente. 

Repetindo o processo muitas vezes (e com o com- 
putador essa repetição não tem significado), acha-se 
grande número de soluções que permitem deduzir a lei 
de distribuição de probabilidade do critério de explorabi- 
lidade (nomeadamente, os seus valores médio e desvio 
padrão). 

Vejamos o exemplo concreto da determinação do 
lucro de uma exploração mineira (definida pelo seu 
valor esperado e a sua variância) à custa do conheci- 
mento do preço de venda unitário (ambos expressos 
pelos valores médios e desvios padrões). 

Cada amostra corresponde à extracção de um nú- 
mero aleatório (entre O e 1), o qual é comparado com as 
distribuições de probabilidade acumulada já conhecidas 
(uma para o preço dae venda e outra para o custo de 
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Preço de venda por ton. 
m= 13800 s= 1583 


Custo de produção por ton. 
m = 9800 s=i%22 


Lucro unitário 


| Amos- | 
tra tt) (2) (3) (4) (5) v 
Desvio normal m+sxXtt) Desvio normal m+sX(3)| L=(2)—(4) 
aleatório aleatório 

0,487 13,89 0,391 | 9,40 4,49 20,16 | 

0,393 12,28 0,939 10,40 2,14 | 4,57 

1,831 16,35 —Q,844 7,97 8,38 | 70,22 

0,243 13,45 —Tt, 43 8,09 | 5,36 | 28,72 
TOTAIS 196,72 | 1038,26 


produção), as quais são do tipo indicado na Fig, 2. 
Assim obtêm-se 50 resultados diferentes para o 
lucro, cujo valor médio ou esperado é: 
196,72 
50 


= 3893/ton 


a o desvio padrão do lucro será: 


f do EUÊ 
/ nt — (SLi)? 


g = V E 


/ 50 X 1038,26 — 196,72º 
eee | 
N 50 X 49 


O intervalo de 90% de confiança é, 
+ 1,64 X 2832 = + 3881. 


antão, 


RR EE 
«4 
3 | 
3 08 | 
E | 
| 
2. | 
E 0.4 
ed | 
D | 
E 0.2 | 
Er ] 
[a 


-3s-25-s O s 25 E 


DESVIOS EM RELAÇÃO AO VALOR MÉDIO 


Fig. 2 


Quanto malor fosse o número de simulações, mais 
correcto seria o resultado, pelo que o uso da simulação 
por computador se torna indispensável em problemas 
complexos, 
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No caso geral de um empreendimento mineiro, o 
método de Monte Carlo estaria apto a tratar problemas 
reais que apresentem grande número de variáveis alea- 
tórias, tal como se refere na Fig. 3. 

O problema é então resolvido à custa da simulação 
de grande número de situações onde, em cada uma de- 
las, é calculado o critério de explorabilidade (por exem- 
plo, o valor actual liquido, ou a taxa de desconto in- 
terna). 


Seguidamente procede-se à determinação da média 
e do desvio padrão correspondentes a esse critério e 
apresentam-se os resultados sob a forma indicada na 
Fig. 4. 

Nestas condições, podem ser equacionadas várias 
alternativas para o desenrolar do projecto, escolhen- 
do-se depois a melhor solução, à custa da comparação 
da distribuição probabilística do critério económico. 

Este critério emprega, pois, a noção de risco, ex- 
presso em termos da variância, maior ou menor, que se 
obtém nos resultados finais. Existirão, assim, dois fac- 
tores para a escolha da melhor solução: 


— o maior lucro previsto; 


— O menor risco, ou seja menor variância da 
grandeza que exprime o critério, 


É costume traduzir esse risco pelo quociente entre 
o limite inferior de confiança e o valor médio, para deter- 
minada rentabilidade, visto assim se avaliar a menor ren- 
tabilidade que é garantida com determinado grau de con- 
fiança. Por exemplo, quando essa probabilidade é 95%, 
o risco será calculado pela relação: 


m — 2e 


m 


Em certos casos, a escolha é simples, como por 
exemplo: 
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PROBABILIDADE 
DE 
OCORRÊNCIA 


CUSTO DO CAPITAL CUSTO DE PRODUÇÃO VENDAS 


TEOR DO MINÉRIO RESERVAS 
Fig. 3 
4 — CONCLUSÕES 
A B 

| Taxa de desconto interna esperada | 25% 20% O estudo econômico de um empreendimento mi- 
| neiro é um processo que envolve três tipos de análise: 

Limite inferior de confiança a 95% 18% | 15% = Ifommácãa: 

— avaliação; 
O caso A é preferivel por satisfazer aos dois cri- — decisão. 


térios de melhor rentabilidade e menor risco. 
Mas pode acontecer que: 
A B 
Taxa de desconto interna esperada 21% 26% 


Limite inferior de confiança a 95% 15% 10% 


PROBABILIDADE 
DE 
OCORRÊNCIA 


VALOR ACTUAL LIQUIDO 


Em geral tal estudo é caracterizado por uma infor- 
mação limitada, que desde início resulta do carácter 
incompleto dos dados de prospecção e pesquisa. Em 
regra, as estimativas de reservas e teores são vagas, 
facto que se traduz por grande incerteza ao aplicar os 
seus valores a cálculos econômicos neles baseados. 

Além disso, a estimativa de custos de produção é, 
na maior parte dos casos, muito dificil de elaborar, re- 
querendo um grupo de técnicos com diferente especiali- 


PROBABILIDADE 
DE 
OCORRÊNCIA 


TAXA DE DESCONTO INTERNA 


Fig. 4 


A decisão neste caso já é complexa e varia de 
empresa para empresa, pois depende da sua estratégia 
de crescimento, de obtenção de lucros, e doutros objec- 
tivos, 

Geralmente cada companhia estabelece essa estra- 
tégia à custa de curvas de utilidade próprias, onde se 
equacionam o lucro (em abcissas) com o risco (em or- 
denadas), procurando um equilibrio entre ambos. 

É pois uma matéria de julgamento que inclui fac- 
tores não económicos, dependendo das estratégias de 
decisão administrativa das empresas, 
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zação e considerável experiência, que se pronunciam 
sobre cada caso concreto. Os financiamentos e proces- 
sos de aquisição de capital introduzem mais dados novos 
que precisam de ser devidamente ponderados e, poste- 
riormente, quantificados. 

Outros factores influentes são também dificeis de 
quantificar: a localização do empreendimento e o custo 
dos transportes, as disponibilidades e o custo da mão-de- 
-obra, a situação política a evolução dos mercados (no- 
meadamente do preço de venda do produto a explo- 
rar), etc, 
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Existe, pois, uma vasta quantidade de informação, 
que corresponde à multiplicidade de hipóteses possiveis 
para o desenvolvimento de uma mina. Para tratar tão 
larga porção de dados e os processar de acordo com 
critérios económicos de rentabilidade, é necessário, cada 
vez mais, recorrer a computadores. Mas tal não basta, 
pois é indispensável tirar partido dos computadores para 
tornar as análises mais perfeitas, sabendo interpretar 
convenientemente os resultados obtidos. 

O critério de rentabilidade dos empreendimentos é 
designado com mais propriedade por rentabilidade espe- 
rada, não podendo dissociar-se da noção de variância 
associada à rentabilidade. 

Há pois a necessidade constante de quantificar o 
risco dos empreendimentos dado que ele faz parte inte- 
grante da noção de rentabilidade. 


INFORMAÇÃO 


AVALIAÇÃO 


| Estimativas dos valores esperados 
dos parametros mais importantes 


| Estimativas da variancia dos pará 


metros mais importantes Monte Carlo 


Imprevisíveis 


Valor actual Tiquido ou | 
| Taxa de desconto interna 


| Simulação pelo metodo de 


Existe ainda um conjunto de dados imprevisíveis, 
que devem também ter lugar dentro do processo de 
estudo da explorabilidade. Exemplos de avalanches, en- 
xurradas, incêndios, nacionalizações e outras imprevi- 
síveis da indústria mineira, podem fornecer elementos 
para julgamento que influenciam também a decisão 
final. 

Esquematicamente, podem indicar-se as fases prin- 
cipais do estudo de um empreendimento mineiro, de 
acordo com Mackenzie, no diagrama que consta da 
Fig. 5. 

Como se observa, a análises completa da explora- 
bilidade de um empreendimento mineiro pressupõe uma 
atitude mental apropriada à indole interdisciplinar do 
problema. 


DECISÃO 


Criterios de rentabilidade 
esperada 


Decisão sobre a ex- | 
plorabilidade do em 
preendimento 


Criterios de variancia 


Fio. 5 
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«1 dos 2 Guindaste eléctricos de 100 tx34m/42tx57m/ 
16tx61m fornecidos para o Estaleiro de Construção Naval de 
Uddevalla, na Suécia.» 


PONTES ROLANTES, GUINDASTES E 
APAR. DE ELEVAÇÃO ESPECIAIS Projecto e fabrico 
TURBINAS HIDRÁULICAS — Fabrico segundo licença de A.C. M. de Vevey, S.A. 
TURBINAS A VAPOR ——————— Fabrico segundo licença de Brown Boveri, Cie. 
CALDEIRAS A VAPOR —————————— Projecto e fabrico segundo licença de Foster 
Wheeler, Co. 


EQUIPAMENTOS E INSTALAÇÕES 
INDUSTRIAIS 


CONSTRUÇÕES METALOMECANICAS MAGUE :... 


ALVERCA DO RIBATEJO = PORTUGAL 
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EXPORTAR 
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Num esforço para aumentar o volume das suas 
exportações, após a encomenda de 200 caixas de 
aço inoxidável para os USA, 


a SOREFAME obteve mais uma encomenda no 
sector de material ferroviário, de valor superior a 
200.000 contos. 


Entre a FEPASA (Brasil) e um consórcio constitui- 


do por 


SOREFAME - Sociedades Reunidas de Fabricações 


Metálicas, S.A.R.L. 


ACEC — Ateliers de Constructions Electriques 


de Charleroi, SA 


MAFERSA Material Ferroviário, S.A. 
INDUSTRIAS VILARES, S.A. 


foi assinado um contrato para o fornecimento de 
30 Unidades Quádruplas Eléctricas. 


A participação da SOREFAME tem especial signi- 
ficado por incluir a execução integral do projecto 
das Unidades Quádruplas, para além do fabrico, 
nas suas oficinas, de 40 veículos motores. 


Os restantes 20 veículos motores e as 60 carruagens 
reboques destas Unidades Quádruplas serão fabri- 
cadas no Brasil. Os motores eléctricos serão forne- 
cidos pela firma ACEC (Bélgica). 


Esta encomenda, obtida em concurso internacio- 
nal, exigiu um financiamento à participação da 
SOREFAME, que foi assegurado pelo Banco de 
Fomento Nacional. 


SOREFAME 
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METODOLOGIA PARA O ESTUDO SÍSMICO DE 
EDIFÍCIOS PREFABRICADOS POR PAINÉIS () 


RESUMO 


A homologação de edificios prefabricados por pai- 
néis em zonas siísmicas requer a realização de vários 
tipos de estudos conducentes à caracterização do com- 
portamento das estruturas dos edificios sob acções 
sísmicas. Tais estudos dizem respeito principal- 
mente a: il idealização do comportamento estrutural 
sob a acção de forças horizontais; ii) análise dinâmica 
em regime linear e não-linear; ilil resistência e ducti- 
lidade dos elementos estruturais isolados ou da estru- 
tura global. No tratamento desses problemas torna-se 
indispensável conjugar de forma adequada a utilização 
de técnicas tanto analíticas como experimentais. 

Nos últimos dez anos têm sido realizados pelo 
Laboratório Nacional de Engenharia Civil (LNEC) diver- 
sos estudos sobre o comportamento sismico de edificios 
prefabricados por painéis, a construir em Portugal e no 
estrangeiro. Sintetizam-se no presente trabalho os prin- 
cipais aspectos da experiência adquirida neste dominio 
e a metodologia que tem sido seguida. 


1 — INTRODUÇÃO 


Nos últimos vinte anos tem-se generalizado em 
muitos países a prefabricação de edifícios por painéis 
de grandes dimensões, o que permitiu reduções consi- 
deráveis no custo e nos prazos de execução dos edifí- 
cios, particularmente significativas em programas de 
construção de habitações de grande volume, 

Adoptam-se presentemente muitos tipos diferentes 
de prefabricação por painéis [1]. Em geral os painéis 
das paredes são feitos de betão simples ou armado, 
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SYNOPSIS 


The agreement for the use of prefabricated panel 
buildings in seismic areas requires different kinds of 
studies in order to determine the major aspects of the 
seismic behaviour of the buildings. Such studies should 
concern: il idealization of structural behaviour under 
horizontal loads; ii) linear and non-linear seismic re- 
sponse; and iii) resistance and ductility of individual 
members and the whole structure. Suitable combination 
of analytical and experimental techniques should be 
used for studying these problems, 

In the last ten years the Laboratório Nacional de 
Engenharia Civil (LNEC) has studied several prefabri- 
cated panel structures to be built both in Portugal and 
abroad. This paper briefly reports the experience of the 
LNEC in this field and describes the methodology that 
is being used, 


incorporando nalguns casos alvenaria de tijolo ou agre- 
gados leves; as lajes são de betão armado ou pré-esfor- 
cado, por vezes com tijolos ou blocos; os painéis e as 
lajos são ligados por juntas de vários tipos. Às alturas 
dos edifícios vão em geral até aos 20 andares. 

Os sistemas de prefabricação por painéis utili- 
zados ou a utilizar em Portugal num futuro próximo não 
foram concebidos inicialmente para resistirem a solici- 
tações sísmicas, de modo que tem sido necessário rea- 
lizar estudos especiais para definir as exigências a que 


(1) Tradução da Comunicação «Methodology for Seismic Studies of Prefabricated Panel Buildings», apresentada ao 4.º Sim- 
pósio Europeu de Engenharia Sísmica, que se realizou no Imperial Collage, em Londres, em Setembro de 1972. 
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estas estruturas devem satisfazer nas diferentes zonas 
sismicas do Pais. O Laboratório Nacional de Engenharia 
Civil tam-se ocupado desta tarefa assim como de emi- 
tir os respectivos documentos de homologação [2 a 5]. 
O LNEC estudou recentemente o comportamento sísmico 
de edifícios de 15 andares, prefabricados por painéis, 
na Venezuela (Fig. 1) [6]. 
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Apresentam-se a seguir, resumidamente, os métodos 
utilizados nestes estudos e as principais conclusões 
obtidas. 


2 — IDEALIZAÇÃO DO COMPORTAMENTO ESTRUTU- 
RAL 


Admite-se correntemente que os edifícios de painéis 
se comportam como consolas verticais complexas, sob 
a acção de forças horizontais [7]. As translações laterais 
totais são devidos às deformações por flexão e corte 
das consolas e à deformabilidade das fundações (prin- 
cipalmente rotação, visto a translação ser desprezivel 
na maior parte dos casos), conforme se indica na fig. 
2 (8). 
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Fig. 2 — Deformada de uma consola sob a acção de forças 


horizontais 
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A importância relativa dos três tipos de deforma- 
ção considerados depende das características dos edi- 
fícios e das funções, sendo em regra necessário proce- 
der a estudos experimentais e a estudos analíticos para 
esclarecer as características do comportamento estrutural 
destes edifícios. 

Assim a realização de ensaios dinâmicos dos pro- 
tótipos sob a acção de vibrações ambiente e de vibra- 
ções forçadas constitui um meio simples e prático de 
determinar as suas deformabilidades e características 
dinâmicas (frequências próprias, modos de vibração e 
amortecimento). 

Os ensaios de vibrações forçadas permitem em 
aspecial determinar a rigidez de rotação das fundações 
[9], conforme se indica esquematicamente na fig. 3. 
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Fig. 3 — Ensaio de vibrações forçadas para determinar a 


rigidez de rotação das fundações 


Com este fim utilizou-se no ensaio de edificios de 
estrutura-parede de 44 andares em Parque Central, Cara- 
cas, um pêndulo mecânico muito simples, constituído por 
uma massa metálica de 2,5 tf suspensa de um cavalete 
(fig. 4). O periodo de oscilação do pêndulo regulava-se 
facilmente, fazendo variar o comprimento dos esti- 
cadores. 
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Fig. 4 — Pêndulo mecânico utilizado nos ensaios dos edifícios 


de Parque Central, em Caracas 
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Por outro lado, os ensaios de vibrações ambiente 
são suficientes para determinar as frequências próprias 
e os modos de vibração, Conhecendo a distribuição das 
massas ao longo das consolas, é fácil calcular a rigidez 
global da estrutura [10 a 15]. A diferença entre a defor- 
mabilidade global e a deformabilidade devida à rotação 
das fundações é igual à deformabilidade devida à 
flexão ao corte nas consolas, consideradas encastradas 
na base, 


É necessário comparar as deformabilidades deter- 
minadas experimentalmente com as calculadas. Um 
problema dificil no cálculo destas deformabilidades é 
avaliar a contribuição dos painéis transversais à direc- 
ção considerada para a rigidez global da estrutura. Con- 
tudo, uma combinação judiciosa de resultados analíticos 
e experimentais permite obter informações precisas acer- 
ca da deformabilidade sob a acção de forças horizontais, 
em regime linear. 

Como se indicou, as frequências próprias e os 
modos obtêm-se directamente dos ensaios de vibração. 
Os seus valores devem também ser comparados com os 
resultados analíticos. 


Em edifícios pequenos e médios (até cerca ce 6 
pisos), a deformabilidade das consolas é muito pequena 
em relação à rotação das fundações. Por esse motivo 
a frequência fundamental e o primeiro modo de vibração 
não se afastam muito dos obtidos considerando um 
bloco rígido assente num solo elástico. Supondo a con- 
sola rigidamente encastrada, as frequências próprias são 
bastante superiores aos valores determinados experimen- 
talmente. 


Os resultados obtidos em estudos analíticos e ex- 
perimentais de uma gama de edifícios de 6 a 15 anda- 
res prefabricados por painéis [9] sugeriram uma fór- 
mula empírica simples para o periodo fundamental, T: 


T=k (1) 
B 


em que H é a altura e B a largura do edifício. Para o 
sistema de prefabricação que foi objecto do estudo 
obteve-se um valor de k de cerca de 0,20 segundos, 
para módulos de elasticidade de compressão não uni- 
forme dos solos com valores compreendidos entre 2,5 
a 5 kgf/cm3. 

A equação (1) tem-se mostrado bastante precisa 
para diferentes tipos de estrutura-parede: painéis pre- 
fabricados, pórticos preenchidos com alvenaria de tijolo, 
e até paredes de betão armado de grande altura, desig- 
nadamente nos edifícios de apartamentos de 44 an- 
dares já referidos. Neste último caso os períodos de 
oscilação fundamental correspondiam a um valor de 
k = 2,25 segundos [11], 

Na fig. 5 apresentam-se as frequências e os modos 
de vibração fundamentais do edifício de 15 andares, 
prefabricados por painéis, a que se refere a Fig. 1 [6] 

A contribuição dos três tipos de deformabilidade 
mencionados é muitas vezes tal que a configuração do 
primeiro modo de vibração resulta aproximadamente 
triangular. 
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DO EDIFÍCIO 
Fig. 5 — Frequências e modos de vibração fundamentais em 


edificios de painéis com 15 pisos [6] 


A grande rigidez dos edifícios de painéis leva a 
que para casos correntes a frequência fundamental seja 
superior a 2 Hz, apresentando-se as frequências corres- 
pondentes ao segundo modo de oscilação bem desta- 
cadas da frequência fundamental e com valores entre 
3 a 5 vezes superiores. Assim, os modos de vibração de 
ordem superior contribuem muito pouco para a resposta 
sismica, pelo que a análise dinâmica pode ser aplicada 
só ao modo fundamental. 

Os ensaios de vibração também permitem determi- 
nar o amortecimento estrutural, que se pode considerar 
do tipo viscoso, O amortecimento relativo para os modos 
significativos oscila entre 0,02 e 0,05, aumentando com 
a amplitude de vibração. 

Finalmente é necessário ter em conta que os pe- 
riodos e os amortecimentos obtidos em ensaios de vibra- 
ção de pequena amplitude são mais pequenos que os 
valores obtidos com vibrações de grande amplitude. 
Como indicação aproximada pode-se dizer que, para a 
análise dinâmica sob acção de sismos intensos, os pe- 
riodos medidos devem ser aumentados de 10 a 20% 
[12]. Atribue-se em geral um valor de 0,05 ao amorte- 
cimento de estruturas de betão armado. 


3 — ANÁLISE DINÂMICA EM REGIME LINEAR E NÃO- 
“LINEAR 


Na maior parte dos estudos sísmicos realizados 
no LNEC as solicitações são idealizadas sob a forma 
de amostras de vibrações aleatórias estacionárias, com 
uma função de densidade de potência de aceleração 
prescrita e uma duração de 30 segundos [13]. A forma 
da função da densidade espectral depende da sismi- 
cidade local e das características do solo. Adopta-se 
frequentemente, para valor médio da densidade na gama 
de frequências O a 5Hz, em solos médios, o de 700 
gal2/Hz. 

Os sismos com intensidade (grau IX na escala 
de Mercalli) provocam esforços muito grandes nas es- 
truturas de painéis, devido à sua grande rigidez. Se a fre- 
quência fundamental for superior a 12 Hz e o comporta- 
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mento for linear os esforços transversos na base atingem 
30% do peso do edifício ou mesmo mais, com uma distri- 
buição das forças sismicas aproximadamente triangular 
com valores máximos no topo do edifício. Com efeito, 
desprezando os moldes de ordem superior, a amplitude 
do deslocamento horizontal X; no piso i é dada por: 

k. A Z, 12) 


em que (7, é a coordenada do primeiro modo de vibra- 
ção, com o valor máximo igual à unidade. 

Supondo que a densidade do espectro de potência 
de aceleração é aproximadamente uniforme no dominio 
de frequências considerado, o valor de À em em para 


para frequências f = 2 Hz é dado por [14]: 


KR, S, |? E — 1 q 
Alc- | 4 | TI É (3) 
or M, 700 0,05 
am que: 
R/M, -——- factor de participação correspondente ao 
primeiro modo; 
S, -— densidade do espectro de potência de 


aceleração (gal Hz !); 


+ — amortecimento relativo. 


Supondo que as forças sísmicas F. são proporcio- 
nais ao primeiro modo de vibração, vem: 


em que « é uma constante e W, É a carga permanente 
ao nível do pavimento i; para que o esforço transverso 
total na base seja C E W, é necessário que se verifique 
Ita: 


& = = [5] 


sendo o coeficiente sísmico, É, dado por: 


a E + (6) 
q EW, g EW, 
em que a,., é a aceleração máxima no topo do edifício 
Supondo que o primeiro modo de oscilação tem 
uma configuração aproximadamente triangular, que o 
peso do edifício está uniformemente distribuído, que 
S, = 700 gal” Hz"! a « = 0,05, da equação (3) 
obtém-se: 


Mis 1 (7) 


e a aceleração maxima no topo tem o valor: 


an=4mt A=0,6 9 (8) 


max 


O coeficiente sísmico fornecido pela equação (6) 
é então: 


C=0,6 20,=0,6xX0,5=0,3 (9) 


Embora os ensaios estáticos e dinâmicos de edi- 
fícios de painéis tenham mostrado a grande resistência 
destes a forças horizontais desde que as juntas sejam 
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devidamente projectadas, as acelerações e as forças 
acima indicadas podem muitas vezes exceder o limite 
elástico, produzindo plastificação em maior ou menor 
escala dos painés e das juntas. 


Nalguns casos é possivel e aconselhável reduzir a 
rigidez do edifício, e portanto a sua frequência própria 
de oscilação, por meio de um primeiro piso flexivel. Este 
método tem sido utilizado em Portugal para edificios 
prefabricados com mais de dez andares (fig. 6). 


Fig. 6 — Edifício de painéis de 12 pisos com andar inferior 


flexivel, Lisboa 


À cedência tende a concentrar-se no piso flexivel, 
protegendo os painéis acima deste. Note-se contudo que 
o piso flexível deve ser devidamente projectado, tendo 
em consideração que nele se concentram os desloca- 
mentos horizontais, devido à descontinuidade na rigi- 
dez. À observação de edifícios deste tipo depois do sis- 
mo de 1969 em Lisboa mostrou que o seu comporta- 
mento estrutural estava de acordo com as previsões. 


Muitas vezes, contudo, é necessário considerar o 
comportamento não-linear dos painéis e das juntas sob 
acção de sismos de grande intensidade. É então neces- 
sário caracterizar o comportamento estrutural sob acção 
de solicitações repetidas de grande amplitude. A técnica 
adoptada no LNEC consiste em ensaiar modelos de gran- 
des dimensões ou partes de protótipos, em vez de pai- 
néis ou juntas individuais. Na fig. 7 indicam-se as di- 
mensões gerais de um edifício de dois andares construído 
por painéis, que foi ensaiado sob acção de ciclos de 
forças estáticas horizontais aplicadas no topo. O edificio, 
apoiado em duas vigas de fundação longitudinais, foi 
submetido a um pré-esforço vertical, a fim de reproduzir 
o peso dos 13 andares superiores. Na fig. 8 apresenta-se 
a fendilhação obtida para dois valores de força horizon- 
tal aplicada. 
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Fio. 7 — Dimensões gerais de um edifício de painéis com 
2 pisos, ensaiado à rotura [6] 


FENDILHAÇÃO PARA F. 428 LF 


Fig. 8 — Fendilhação em painéis longitudinais 


Os diagramas forças-flechas obtidos nestes ensaios 
foram adoptados na análise dinámica não-linear, utili- 
zando o método da consola equivalente [16]. Este mé- 
todo consiste essencialmente em idealizar consolas com 
propriedades dinâmicas idênticas às da estrutura real. 
A sua deformabilidade por flexão e por corte pode ser 
representada por diagramas momento-curvatura B es- 


MT, b 
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M - Momento flector 

T - Esforço transverso 
B(M)-Curvatura 
P(T)-Rotação 


bilingar com 
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transições parabólicas 
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forço transverso-rotação de forma bilinear com transi- 
ções parabólicas (fig, 9). A análise dinâmica das con- 
solas equivalentes baseia-se numa integração numérica 
passo a passo, pelos métodos de Runge-Kutta [17]. 


O método da consola equivalente é muito geral e 
pode ser aplicado a estruturas de pórticos e estruturas 
mistas. Recentemente foi utilizado para optimizar o pro- 
jacto dum andar flexivel num edifício de painéis, de 13 
pisos, a fim de manter dentro de limites adequados as 
forças sísmicas devidas a sismos de grande intensidade, 
Foi assim possível reduzir a 0,15 o coeficiente sis- 
mico [18]. 


4 — RESISTÊNCIA E DUCTILIDADE DE ELEMENTOS 
ISOLADOS E DA ESTRUTURA GLOBAL 


A influência do comportamento não-linear na res- 
posta sísmica traduz-se correntemente por factores de 
ductilidade que medem a redução das forças sísmicas 
ou dos coeficientes sísmicos provocada pelo comporta- 
mento não-linear, A ductilidade dos elementos indivi- 
duais é sempre superior à ductilidade da estrutura no 
seu conjunto, dependendo ambas da intensidade sis- 
mica e do tipo de estrutura. O problema tem sido dis- 
cutido por muitos autores, admitindo-se duma maneira 
geral que os factores de ductilidade global podem ter 
valores até 3 em estruturas de pórticos, em que se 
poda formar um grande número de rótulas plásticas nas 
zonas sujeitas a maiores esforços. Em estruturas de pai- 
néis, em que a cedência é devida principalmente a fen- 
dilhação por tracção e por corte nos painéis e nas 
juntas, os factores de ductilidade são mais pequenos 
que nos pórticos. 


A forma mais directa e precisa de determinar a 
ductilidade global consiste em executar análises dinã- 
micas lineares e não-lineares e comparar as forças sis- 
micas obtidas. Este método pode muitas vezes simpli- 
ficar-se, atendendo ao facto bem conhecido de que, 
para valores baixos de ductilidade dos elementos, os 
deslocamentos sísmicos não são muito afectados pelo 
comportamento não-linear, Podem-se assim calcular os 
deslocamentos máximos por meio de análises dinâmicas 
lineares e comparar os deslocamentos obtidos com as 
deformadas devidas aos diferentes valores das forças 
sismicas em que se teve em conta o comportamento 
não-linear. Esta técnica foi aplicada no estudo do edi- 
fício de painéis de quinze andares já referido [6]. 
Apresentam-se na Fig. 10 os resultados obtidos, sob a 
forma de deformadas para distribuições triangulares das 
forças horizontais, para comportamentos linear e não- 
“linear, correspondentes a diferentes valores do coe- 
ficiente sísmico C. O comportamento não-linear iniciou- 
-sa para C = 10,10. As deformadas correspondentes a 
C > 0,10 foram obtidas por combinação dos resultados 
dos ensaios de rotura, mencionados no parágrafo ante- 
rior, com a deformabilidade elástica determinada em 
ensaios de vibração ambiente. 

A intensidade sísmica máxima considerada provo- 
cou deslocamentos de 2,4 em no topo do edifício, que 
correspondem a um coeficiente sísmico da 0,32 consi- 
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Fig. 10 — Deformadas de um edifício de painéis de 15 pisos 


considerando comportamentos linear e não-linear 


derando um comportamento linear. Ássim a ductilidade 
global neste caso tem o valor 0,32/0,25 = 1,28. 

Está actualmente em curso no LNEC o estudo do 
comportamento sísmico de edifícios de painéis, de 8 
a 20 pisos. Os resultados já obtidos confirmam que a 
ductilidade alobal, sob a acção de sismos fortes com 
as intensidades indicadas no parágrafo 3, oscilam entre 
1,2 0 1,5. 

A resistência às forças horizontais e a ductilidade 
última dos edifícios prefabricados por painéis são muito 
influenciadas pelo tipo dos painés e pela pormenori- 
zação das juntas. Assim, cada tipo de edifício exige 
normalmente estudos experimentais específicos a fim 
de determinar essas caracteristicas. 


5 — CONCLUSÕES 


A experiência obida nos estudos sismicos de edi- 
ficios prefabricados por painéis mostra que este tipo 
de estrutura pode ser utilizado em condições econômicas, 
mesmo em regiões de forte sismicidade. Contudo, cada 
novo tipo de construção exige a realização de estudos 
que definam as condições que os edifícios devem satis- 
fazer em função da sismicidade local, Estes estudos 
— tanto experimentais como analíticos — devem satis- 
fazer aos princípios da engenharia sismica moderna, 
permitindo comparar os deslocamentos e esforços sis- 
micos que se desenvolvem na estrutura com a respec- 
tiva ductilidade e resistência. É desta comparação, cujos 
termos devem ser quantificados sob o ponto de vista 
estatístico, que resulta em última análise o julgamento 
global sobre a segurança da estrutura relativamente às 
solicitações sismicas, 
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SOLUÇÃO DA EQUAÇÃO DE POISSON 
PELO MÉTODO DE RITZ 
- PROCESSO GERAL PARA QUALQUER CONTORNO 


RESUMO 


E aplicado à equação de Poisson bidimensional com 
condições na fronteira de Dirichlet a técnica da semi- 
-discretização. Mostra que o sistema de equações dife- 
renciais obtido permite uma formulação variacional 
equivalente. A solução do problema variacional pelo 
método de Ritz, reduz-se ao cálculo de um conjunto de 
funções que se exprimem explicitamente em termos de 
coordenadas de pontos sobre o contorno que limita o 
dominio da solução, 


INTRODUÇÃO 


Os métodos numerados de integração da equações 
tipo elíptico baseado no cálculo das variações, têm 
sido usados com sucesso. Desses métodos os designa- 
dos por semi-analíticos (de Ritz, de Kantorovich, etc.) 
possuído a vantagem de uma formulação geral, apre- 
sentam inconvenientes na aplicação prática sempre que 
o contorno do domínio de solução do problema não 
tem uma descrição analítica simples. 

No que se refere ao método de Ritz, apesar dos 
resultados obtidos na sua extensão a contornos de ex- 
pressão analítica complexa,* não se pode considerar, 
do ponto de vista prático, de aplicação geral, uma vez 
que sempre que a forma do contorno não é simples, 
a forma da solução proposta vem muito complexa e o 
cálculo das constantes torna-se excessivamente tra- 
balhoso. 

Os inconvenientes apontados podem eliminar-se, em 
vez de se procurar a solução da equação às derivadas 
parciais, se determinar a solução do sistema de equa- 
ções diferenciais ordinárias que se obtêm a partir da- 
quela pela aplicação do método da semi-discretização. 

A solução do sistema é a resposta ao problema e a 
sua precisão é dependente do intervalo de discretização 
utilizado. 


* Kantorovich e V. Krilov — «Aproximate Methods of 
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SUMMARY 


The straight line method is applied to bidimensional 
Poisson equation with Dirichet boundary conditions. 
tt shows that the resulting system of differential equa- 
tions allwos an equivalent variationa! statement. The 
solution of the variational problem by Ritz method is 
reduced to the calculation of a set of function which 
are closed forms in terms of the boundary point coordi- 
nates. 


1. EQUAÇÃO DE POISSON 


Considere-se o domínio bidimensional D limitado 
palo contorno + a admita-se r tal que qualquer recta 
paralela ao eixo dos x unicamente o intersecta em dois 
pontos, 


Na direcção y, r é limitado inferiormente pela 
recta y = O a superiormente pela recta y = 1 
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Põe-se o problema da determinação da função 
vix,y) que satisfaz no domínio D a: 


atv ay 
Es ayt glx,y) (1 | 


g às condições de Dirchlet 


vaglxy) para xy em T (2) 
Suponha-se 
gix,o0)=wlx) 
( õ (3) 
e lx, l)= vlx) 


a considere-se a transformação: 


Us velM= vg) É vg (4) 


A equação (1) toma a forma 


& as condições na fronteira (2) vêm: 


uz e (x,y) = (= Mo) = von tP ley) par xy, emr (6) 


As condições impostas no contorno + em y =D 
ay=1 são 


P(x,0)=0 


Pix, 120 ds 


3 — SEMI-DISCRETIZAÇÃO 


Defina-se um intervalo h, pequeno, da variável y 


he (8) 


e considere-se o conjunto de rectas 


yry ih LISO, ba Éava n, n+1) (9) 
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Sejama, eb, (i=0,1,2...nn+1) as abcis- 
sas dos pontos de intersecção das rectas y — y; (i=o, 
1...n+ 1) com 7. 


Designando por u; o valor da função incógnita 
na recta y = y; 


uj=Ulx,y;)= ujlx) (i=0,1%,..nn+ 1) (10) 


com w=0 6 Un +10, ter-se-á: 


a2y 


Do Ea dC o pm 
E 2a h2 2 tgyé y=Y) 
(= 1,2 sm 


Fazendo h suficisntemente pequeno, para que se 
possam desprezar os termos de ordem igual ou superior 
a h' em (11) e pondo 


fi fix, yjl E Ta És 1) (12) 


a equação (5) aplicada em cada uma das rectas y = y, 
(i=1,2... n) toma a forma: 


dêuy -2Uyouz | 
dx? 


2 
| dºu us -2U.+eU, 

E (i=zZ, 3 nas n-1) (13) 
É de | 


d?un Un-1-24n 
e dy rim 


f 
tgx? hº i 


Em escrita matricial é: 


dx? 2 (14) 
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onde 
u2 o 4 d | 2 
O 1-2 1 1 
| a 
Un fn 


As condições aos limites a que u tem de satis- 
fazer são: 


usP-Q(a;,y;) Para qa, 
a ii (Gel, É we B) 


uz 9: P(b;,yj) 420) 


para 
(15) 


3 — FORMULAÇÃO VARIACIONAL 


Como se demonstra no Apêndice A, a solução 
da equação (14) conjuntamente com as condições (15) 
é equivalente à pesquisa do vector u que, satisfazendo 
a (15), minimiza a funcional 
* x 


ot: 3[3] (UU -—; UA us2 Fu) (16) 


onde: à é o operador que actuando sobre um vector r 
qualquer realiza a operação 


— [J] é o operador que actuando sobre um vector 7 
realiza: 


— [U] é a matriz de elementos 


- E. K 6, vio U e Ud u 


4 — METODO DE RITZ 


Considere-se o conjunto de Wp (x) [k=1, 2 
- m) independentes e satisfazendo a: 


Wiilaj)=0 . 
(o Lo É 064% 1) (17) 
Wuilb;)=0 


B com componentes tais que 


Wear, Man e Wed Mi 


sejam de ordem h. 


m 
Exprimindo u por Us Sind CW, (18) 
E 


onde C, (j=1,2 ... m) são constantes e S(x) sa- 


tisfaz a 


promo 


(19) 


a expressão (18) satisfaz automaticamente às condições 
aos limites. 
O valor aproximado da solução u, obtém-se pro- 


curando os c (j= 1,2... m) que conduzem ao menor 
valor de (u), quando se substitui (18) em (16). 
Isto é: 


| a Th m m 
Blevereo- eme SUS EM (8-1 [S + EM] 
m m 
Al(s.Z c, wzF(seLc, wi)) : 


m 
«s[5)] Ss. SAs-2Fs-2L c, (Sms SA wi Fwk) 
"EE, (Wi wWj = Tg Wi A Wi ] (20) 


onda 


S:[:65;] Wir [cidi;] 


Impondo as condições de mínimo a O leyca, Em) 


22.0 (k=1,2...m) (21) 


k 


pat] 


* As matrizes serão representadas por maiúsculas em tipo grande, os vectores por minúsculas em tipo grande. 


** ma expressão encerrada por [ 7) representa uma matriz. Uma expressão encerrada por ( | representa um vector 


... é — Simbolo da Kronecher 
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obtém-se 


UMES; [Wiwj 15 MAM) Mol; S Ao F Wi js 


Ju 
lk= 1,2... m) (22) 


donde, os C, (j= 1,2 ... m) são determinados pelo 
sistema de equações algébricas 


fm 
PA Ly; Cry lie= 1, 2... mi) (23) 
jat 
onda 
n n n n=1 
sy Pi E (Ea Chjiioa "2 À Chii* fi Ohiisio!) 
(24) 
com 
Di 
ai 
bi 
ejist / Wit Wir dx (26) 
Qi 
e 
a. e bi bi 
Ys-ç WE: dx — E Wi dx-2 SW: OX 
E ad ki pElh, d, |=1Pki É Hj [ki 
bi bi 
N=l n 
pod E j.l Wi dx Ff tw “dx | 
a de ) ja J, io (27) 


Estabelecida a forma das funções S, (x) e Wy (x) 
(i=1,2...mek=1,2... m) e conhecidas as fun- 
ções f (x) (i= 1,2... n), o cálculo dos integrais que 
figuram em (25), (26) e (27), podem realizar-se anali- 
ticamente pelo que os coeficientes Qkj O Yk do sistema 
de equações (23) serão expressões conhecidos dos 
valores a, e b, (i=1,2... n). Nestas condições o 
cálculo realiza-se duma forma geral para qualquer 
contorno. 

No Apêndice B, são dados valores de (25), (26) 
e (27) para vectores S e W,. que satisfazem as condi- 
dições gerais impostas. 


5 — EXEMPLO 


O método proposto foi ensaiado no caso particular 
da equação 
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com u= O paraxey satisfazendo a x'+ y = 1. A solu- 
ção exacta 6 u = '/, (1-x*-y;) 

Utilizando um número de intervalos n + 1 = 50,10 
constantes arbitrárias e recorrendo às expressões apre- 
sentadas no Apêndice B, obtiveram-se os seguintes 
resultados: 


[y] x | solução obtida | solução exacta | 


| Vi-x? 


u dy dx=.3927 (valor calculado) 


“ 
bi 


n 
h | uiaxeanis 


(valor exacto) 


O tampo de cálculo consumido no IBM 360/44 
foi de 49 seg. 


6 — CONCLUSÕES 


— O método é de aplicação geral na solução da 
equação de Poisson bidimensional com condições fron- 
teira de Dirichlet. 

— É possivel obter precisão satisfatória com eco- 
nomia de tempo de computação e utilizando computa- 
dores de capacidade de memória limitada. 

— À extensão a outros tipos de equações e a pro- 
blemas tridimensionais parece possível. 
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APÊNDICE A 


Considera-se que 


6 (ue S[s)[Uu-o; UAuez Fu) (A 1) 


tem um valor minimo o vector u que satisfaz as con- 
dições fronteiras. 
O vector 


com dla,)=0 à 


E a pr 
fronteiras. 


n. 1) de ordem h, satisfaz as condições 


Se g (U) é minimo, a 6 (Ucad)s0 para «=0,. 


da 


Substituindo (A2) em (41) vem: 


o lusa di=8 [s](Uu'.za Du-aíDd'-., (UAu-2aDAuU 


. *DAd)-2aDf) (A 3) 


onde D = [d,ô 1. Recorreu-se a UAd = DAu porque A 
é simétrica, 

Derivando (A 3) em ordem a a« e igualando a zero 
a expressão obtida para a = O, vem: 


8[u)[2Du-.2-DAu-2Df).o (A 4) 


O termo 


8 [3][2Du') 


transforma-se por 


n má n a | nd; 
2 q is maio PÁ ue Idjujldx-2 À ) d; uídx= 
db; 


2 x [d; vila, «4 8[4] (D vu? -23 [SJ/D u”? 


Substituindo (A 5) em (A 4) obtém-se: 


sa (D (ue Au-fll-o e 
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Como D é arbitário pode tomar-se 


d,+0 e 


k d=0 (1 nai La 2 .. k-1 ' 


Ke to 


“— n) é fazer dj (x) = O no domínio 


[ax x- A: xÁ, by 


com A pequeno de forma dy (x) toma o mesmo sinal 
em [x-4, x + A] nestas condições (A 6) toma a forma 


x À 
| atcd (Au) -tdáxco 


xd (47) 


o que implica 


E (A 8) 


O que foi feito para a componente k pode fazer-se 
para as restantes componentes e para todos os valores 
de x no domínio considerado, donde 


2 
dl Us; ÁU=f 
dx h (A 9) 
é válida em todo o domínio. 
APÊNDICE B 
Considere-se 
qi- Pi 
sjlx)= LA (x-aj ep; [1= 1,2 eh] (B 1) 
[io 
que satisfaz a 
|6i lajl=p. 
| 1=48 son, 
silb;)=q; (B 2) 
Er a a L=0;' 
Wi! Em (id. sin (KT io) 
que satisfaz a 
w.-la,)s O 
Ms EE Ls É mes Hi) 
We 19,120 (B 3) 
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e tome-se 


Suli) tal que: Oto )=Gpln+I)z0. 


O vector 


m 
JS +53 cjW) 
J=1 (B 4) 


satisfaz as condições indicadas no texto. 

As funções propostas (Bl e (B2) conduzem às se- 
guintes expressões para as funções gerais indicadas 
em 4.: 


bi 2,4: ; Sm O se kéj 
Ls k 1-0; ; ú Gi; ? 2 
Dusio auiLBjlil LL, .cos LT ——. |, cos[ jTL —— oe] 2 TUR ” 
ejrJ i lbi-aj) ( rar) b; a) apt) Fora; para ks) 


(B5) 
Pi 1-8 1 E-a 
Seta) 9 [Hg Mt Domin [ JTE Sa, ) Si (TE o da 
Op UILO it g-tsintpbj-q)-sin (pajeg = golsintebym) =sintropes ih Isa pf) 
ediaorçásii lua pe) (B6) 
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onde 
paTllkibe-ae)=)lby= 0,1) 0b,=mj)/lby-04) 


quTljaçlb;-ai-ko lby-a 1) /lb; -m 14 iby-m;) 


raTElkiby-ogl=jlb; sm; tb;=0])4 lby-a4) 


s=-Mijogib;-alska; by -m q)! Cb; - e)! Uby-m |) 


A q, *P: b;-a| n QrR; hs dj bico; 
Esuls nã q 4 ay =. 2F alii À 
K” pot ER dh By KTL po k KIT h k Dial “Biel k Í 


E-0) : 
aytid fo tysin (NE pet? para k impor 
aj 


Breç)St para 4 pa (87) 


5 
l Gero) [asin[k TE 


Os valores de ep (i) podem tomar-se como 
fjtk) folk) 
O lida - Ene1=i) (B8) 


com f. |k) e f. (k) números reais maiores que zero. 
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Mm ires 
| 632086 
(7 235670. 
WZ see 

Big! 


gpa 


Uma sugestão: Solicitem a nossa 
Assistência. Desta forma não terão 
preocupações. Basta um telefo- 
nema. O mesmo poderão fazer se 
precisarem urgentemente de uma 
peça subressalente. 

A propósito: pelo aumento da pro- 
dução e por medidas de racionali- 
zação foi-nos possível reduzir, con- 
sideravelmente, os nossos 


A construção dos aparelhos diferen- 
ciais eléctricos DEMAG não é uma 
solução de compromisso, o que está 
provado pelos milhares de aparelhos 
em serviço, em mais de 100 países, 
desde há dezenas de anos. No en- 
tanto, mesmo o melhor produto ne- 
cessita de alguma assistência, com 
determinados intervalos. Assim, por 
exemplo, e segundo as normas de 


prevenção de acidentes, os aparelhos DEMAG prazos de entrega. 

deverão ser inspeccionados uma vez por ano. MM 

É evidente que se gasta um ferodo, ou que deve Queiram enviar-nos o talão (em baixo, do lado 
o mesmo ser reajustado. Igualmente, deverão ser ME direito). Teremos muito gosto êm informar- 
verificados, com regularidade, os cabos e as suas W -vos acerca dos nossos serviços de 

guias. Vs. próprios poderão executar estes traba- | assistência ou do nosso 

lhos, ou utilizar os nossos serviços de assistência. "programa de fabrico. 
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IDENTIFICAÇÃO DE SISTEMAS LINEARES POR 
AMOSTRAGEM FINITA DO DOMÍNIO DO TEMPO () 


RESUMO 


Apresenta-se um método para identificação de 
funções caracteristicas de sistemas descritos por egua- 
ções diferenciais ordinárias lineares de coeficientes 
constantes a partir dum número finito de valores da res- 
posta impulsiva tomados em pontos equidistantes no 
dominio do tempo. 

Indica-se a estrutura dum programa para o cálculo 
automático dos polos e zeros da função caracteristica 

a identificar a partir dos valores observados. 

O método é aplicado à determinação experimental 
da função de transferência dum filtro LCR passa-baixo 
de ordem 5. 


1. INTRODUÇÃO 


É conhecido [1] que um sistema cujo modelo ma- 
temático seja uma equação diferencial ordinária linaar 
da coeficientes constantes pode ser caracterizado por 
uma função característica H(s) dada por 


Ns) 


H(s) = D(s) (1) 


em que N(s( e Dis) são polinómios em s de graus 
9n º 9p respectivamente. 
Para sistemas de largura de banda limitada, tem-se 


90 >9n +1 (2) 


em que o valor de gp identifica a ordem do sistema. 

Designando por y(t) a resposta do sistema a uma 
função ult) e por Y(s) e Uls) as respectivas trans- 
formadas de Laplace, representam-se na figura 1 as 
relações entre as diferentes grandezas, quer no domínio 
do tempo, quer no domínio da frequência complexa. 

O objectivo do presente trabalho é a identificação 
da His) pelos seus polos e zeros (raizes de Dis) e 
N(s) respectivamente), a partir duma colacção finita 
de valores y(t) tomados a intervalos iguais, para o caso 


J. M. VIEGAS GONÇALVES 
Eng.º Electrotécnico |. 5. T. 
Assistente do |. 8. T. 


ABSTRACT 


This paper presents a method for characteristic 
function identification in systems which are described 
by sets of linear ordinary differential equations with 
constant coefficients, based on finit sampling in time 
domain of the system's impulse response. 

The setructure of a program used to compute the 
location of poles and zeros from the sampled values is 
described. 

The method is applied to the experimental deter- 
mination of a fifth order LCR low pass filter transfer 
function. 


em que y(t) representa a resposta impulsiva do sistema. 
Para isso estabelece-se a equação de recorrência a que 
vit) deve satisfazer e a partir desta, obtém-se os polos 
a os zeros de His), 


()=h(t)ult) = 
=/h(t-gulydr 


e] 


Yís)= Hs) Uls) 


Fig. 1 — Relações funcionais entre as grandezas 


2. TEORIA 
2.1 — Resposta impulsiva 


Quando H(s) tem a forma indicada em (1) por 
decomposição em fracções parciais, podemos escre- 
ver [1] 


(+) Publicação do Laboratório de Medidas Eléctricas do |. S. T., integrada no Projecto de Investigação TLE-5 subsidiado 


pelo Instituto de Aka Cultura 
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Hgj)= 2 |— 


| 
RE 


(3] 
(S—p,) É 


em que p; (i= 1,2, , m) são as raízes diferentes do 

polinómio D(s), com graus de multiplicidade q; 
m 

( 2 M=99). 

[= 1 


A resposta impulsiva toma então a forma 


Cipi | 
E sr o id E 


vt = htt) = 
i 1 ls (pus 191 


(4) 


Para o caso em que vit) é conhecido através dum 
conjunto de valores igualmente espaçados no tempo 
de 4 t teremos 


m CA tt = 1 
w ' 
yinat) =y, o En 
i=1l 1 
nt; — 1 ef Mt pn (7) 
n=0,1, 


A determinação dos valores de P; e C, poda enca- 
rar-se como um problema de interpolação da função vit) 
por funções do tipo exponencial [2] quando os P, 
são todos diferentes. 


2.2 — Equação de recorrência 


Para os sistemas em causa os valores das gran- 
dezas medidas a intervalos de tempo iguais, satisfazem 
a uma equação de recorrência de ordem 96. 

Com efeito, a solução y(t) pode considerar-se 
como a combinação linear de gy funções — fi (t), 
É (E), cs ; fas (t) — linearmente independentes. 


Definindo as seguintes grandezas: 


Yitdo= prooylt) À = A, 
vo tt) As 

o Ds (D | As 

EF (t)I1 = fi (t)... fan dt) 7 


fo tt)... tão ft) 
EgD— 1) (gDb—1) 
f,9 PES as Ego (t) 


taramos 


Y(ty)y=[F(t)] A, A=T[F(MIl-A YA) (8a, b) 
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Conhecido o valor de Y lt) num certo instante, 
isto é, conhecidos os valores da função e das suas 
tg, — 1) derivadas num certo instante, podemos deter- 
minar as constantes A, por (8 b). [F(t)] não é sin- 
gular, pois por hipótese, os f, são linearmente inde- 
pendentes. 

Para o instante t + à t podemos escrever 


Titratm=[Fi(t+rati] A 
ou ainda, atendendo a (8 b) 


Yit+ram=[Fltrad) CF(dI-Y (= la) Ytt) 
(9) 


em que [a] = [Flt+ at] [Fit)] 1 


é independente de t, pois 


ví 


Fig. 2? — Mudança de origem dos tempos 


como os f, não dependem do instante tomado para 
origem dos tempos, referindo-nos à figura 2 podemos 
ascrever atendendo a (9) 


Ya= [EF t+At] [E (U]-lY, 
ou ainda 

Y= [Fltt+rAtU] CFM] o 1Y, 
a portanto 


[F(tt+rat] CFAOIJ-T=- [Fit+at) CF(t)]— 
= [a (A t)] (10) 


A expressão (9) permite-nos calcular os valores da 
solução e das suas derivada no instante t + At como 
uma combinação linear dos valores tomados no instante 
t, qualquer que seja t. Então, os sucessivos valores da 
solução e das suas derivadas tomadas a intervalos At 
satisfazem a uma equação de recorrência de ordem g,. 
Para o mostrar segue-se o processo indicado em [3]. 

Considera-se a equação característica da matriz 
[a (4 t)] 
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det ([alát)] — à [1]) = O 


Z ajpi=0 (ago = 1) (11) 


Como pelo teorema de Cayley-Hamilton a matriz 


[a (4 t)] 


satisfaz a sua própria equação característica, vem 
90 | 
5 a [a(at)li= [0] 
|=0 


Multiplicando à direita ambos os membros por 
Y [t) e atendendo a (9) obtém-se a equação de recor- 
rência pretendida 


90 
5 aYittristi=o (12) 


|== 0 


Tomando em (12) apenas os valores da solução 
y(t) e designando por y, os valores de ylt) para t 
= nÃt podemos escrever 


O) n=0. os U8 


Para a equação de recorrência (13) a solução geral 
Yn pode escrever-se [1]: 


Plano dos À, 


m 
= 2 (Apto + Ag pe an 
j=1 


n=0, 1, «.. 
(14) 


em que os h; são as raízes de grau de multiplicidaden, 
da equação característica da equação de recorrência, 
a qual toma a forma (11). Os à; representm, portanto, 
os valores próprios da matriz [a (A t)J. 

Identificando (14) com (7) obtém-se 


Ci Att 

Ag = + eg = 15 a) 
; "o tg] 

dy = eMAt (15 b) 


relações estas que nos permitem obter H(s) conhecida 
a solução geral da equação de recorrência. 

Deva notar-se que a relação (15 b) não define 
univocamente o valor de p, dado o valor de X,. Com 
efeito tem-se 


1 2ka 
= log +] -— =. 
P; At i 


ue Ac Ti vai 
àt 


(16) 


em que [Im 4 log Ai k | €< q resultado este que se on- 
contra representado na figura 3. 


Plano dos p, 


1 


Fig. 3 — Transformação p =. log» 
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No caso em que o intervalo de amostragem for tal 
que se verifique 


E q 
At ET =— (17) 
2 “ max 
em que 
asa EMA E tm OM] = ti Ds th 
antão 
Im (Pi) | < — 
At 


a vem sempre k = 0 em (16). 1 

Não se podendo garantir (17) é ainda possivel 
determinar a parte imaginária de p, realizando dois 
ensaios com intervalos de amostragem At, e át; inco- 


mensuráveis e determinando os valores inteiros k, e k- 


(Únicos) para os quais 


a ts. E a us = & qm Jo 
At à t 


(18) 


Para exemplificar o que foi exposto, consideremas 
um sistema de segunda ordem, de que se conhace a 
resposta impulsiva — Figura 4. 


Fig, 4 — Resposta impulsiva dum sistema de 2,º ordem 


Começa-se por determinar os coeficientes da equa- 
cão de recorrência (13) utilizando ,por exemplo, as 4 
primeiras observações 


“doWtrta y+Hy=O 


“Boy T+ a; Y +ty=0 


obtiidos os valores de a, 6 a, determinam-se os valores 
A O A:, raizes de (11) 


ão + a; A+HN=O0 
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A partir de (15 b) podemos já determinar p, e p: 


1 1 


= - Log À = —— LOG À 
Pp: AA g A Pa Rr 9 As 


Para determinar os valores das constante À, e À; 
da solução geral da equação de recorrência, não temos 
mais que utilizar (14) para as duas primeiras observa- 
ções, por exemplo, e, considerando A; E As 


A. + Ay = Vo 
As Mm + As My 
Os valores de C, e C. dados por (15 a), coincidem, 
neste caso (hM > A:), com À, e À; o que nos permite 
escrever 


A, A. | 


Hits) = - + - 
S—pP, S—P, 


podendo ainda passar à forma (1) reduzindo a expressão 
anterior 


SYo TS (As P, + Às P,) 


Mi = DO Refil 


No caso em que y = O vem 


A, (P; — Ps) 


Gt E 


3. APLICAÇÃO 
3.1 — Descrição do programa de cálculo 


Com base na teoria exposta elaborou-se um pro- 
grama para o cálculo automático dos polos e zeros da 
função característica, a partir dos valores de N obser- 
vações afastadas no tempo de At da resposta impulsiva 
dum sistema de grau g, . Supõe-se que At verifica a 
relação (17) e ainda, para que o cálculo seja possivel, 
quo N> 2 go. 

Descrevem-se a seguir as diferentes fases do pro- 
grama, 


| — Determinação da equação de recorrência 


Como as observações de que se dispõe vem afec- 
tadas de erros experimentais em lugar de (13) teremos 


90 


2 A, Yn+i dn 
[= o 


DO da sos NB =" 


Os valores dos a, são determinados de modo a que 


35, seja mínimo [4]. 


1) “Corresponde ao teorema de Nyquist para o caso em que vit] é periódico limitado na frequência a wmax 
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O sistema de equação a que se chega é resolvido 
pelo método de redução de Gauss [5]. 

Se o valor de g, tomado for maior do que o 
correspondente ao do sistema em estudo — 80, — 


matriz do sistema de equações seria singular. Tal não 
acontece quer devido aos erros das observações (3, 0) 
quer aos erros de truncatura nos cálculos, obtendo-se 
um sistema de equações «mal condicionado». 

Para atenuar este inconveniente, os cálculos são 
programados em «precisão dupla» (16 dígitos signifi- 
cativos) nesta fase. 


|| — Determinação dos polos 


Com os valores dos a, obtidos calculam-se os À, 
a partir de (11). Utiliza-se o método de Nickel [6] 
para a determinação das raízes da equação caracteris- 
tica. Obtêm-se os valores dos p;, através de (16) su- 
pondo k = o. 


Hi — Determinação dos zeros 


Com os valores de à, calculados obtêm-se os À, 
aplicando (14) às gp primeiras observações, No caso 
em que 9, > o, os coeficientes A, em número de 
(9p —89,) serão muito menores do que os restantes, 
na medida em que os erros das observações sejam pe- 
quenos em relação aos valores destas. 

No caso em que os X; são todos diferentes os C, 
coincidem com os A,. 

A partir da forma indicada em (3) passa-se à 
forma (1) calculando-se as raízes do polinómio que se 
obtiver para N(s). 


3.2 — Exemplo de aplicação 


O programa anterior foi escrito em Fortran IV e 
ensaiado num computador IBM 360/44 do C.C.U.T.L,, 
utilizando exemplos numéricos e experimentais, Descre- 
ve-se a seguir um dos exemplos experimentais, bem 
como os resultados obtidos. 

Projectou-se um filtro Tchebyshef passa-baixo cujo 
esquema de encontra na figura 5 e onde se indicam os 
valores medidos dos diferentes elementos utilizados na 
sua realização ('), 

Nesta caso, teremos 


U, (s) =H (5) Us (s] 
em que 
A 
H (8) = 


Para obter u, (t) usou-se um gerador de onda 
quadrada com frequência suficientementa baixa, para 


S—P)(S=—P) (S5—P (8=?Pj) |8=P 


800 9 195ômH 250  1958mH I5A 
UOL “000” Ea 
utt) esa |u(t) 
| 09nF 153 nF WanF = 
Fig, 5 — Filtro- Tchebyshef N=5 1 db fç = 5209 Mz 
Ro = 600 0 


que cada meio periodo de wu, (t) possa ter interpretado 
como a resposta à função escalão — u,, lt) — (fig. 6). 


uft) 


Fig. 6 — Resposta à onda quadrada e resposta à função escalão 


Então, 


Us, (s) = His) E 


Aplicando a teoria exposta utilizando u,, (t) — res- 
posta à função escalão — obtemos os polos e zeros de 
H(sl/s. Aparecerá a mais um polo na origem 
(P gp+1 = 0), relacionando-se os coeficientes obtidos 


para os polos com os coeficientes correspondentes a 
His) por: 


C, [H(s)] 
C, [H(s)/s] = De = Ta Es sm dg (19) 
E 
| No) 
Con + [H(s)/s] = % 
" (-p;) 


Assim, como no exemplo em estudo H(s) é de 
grau 5 a análise foi feita com grau 6. 

Para obter os valores de u, (t) em intervalos es- 
paçados de At utilizou-se uma unidade de amostra- 
gem (*) segundo o esquema da fig. 7. 


“ O valor referente à resistência da bobina correspondea F = 5 kHz 


(.) «Sample and hold circuit», 
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Filtro 
em 
estudo 


Comando | Unidade | 4281) 
de de | 


Voltímetro q.m. 


Y-Registador XY 


atraso | amostrager 


Contador 


alectrônico X-Registudor XY 


Fig. 7 — Esquema para a determinação de valores de u, (t) em 
intervalos espaçados de At 


O ensaio foi feito com uma frequência do gerador 
de 300 Hz tendo-se registado valores com At = 25 us 
O tempo de amostragem era de 3 us. Na figura 8 apre- 
senta-se um gráfico tal como foi obtido com o regis- 
tador X Y. 


E ] 7 [alça i 
É F A as Rr I 

- tu prdir Prato pia a 

! 1 IE | À AN | | 


 FILIRO TEMERYSMEF  NaS tdh Easponda o 
RESPOSTA À PUNÇÃO EStALÃO 


| | E! “H 


| 
| 


| | 
| 

a 

; q hide Ir | 

] | o há BE rm pEd da E dia 

| | hi -CHESHES Era at tio, a] ami 

; : ! | 

1 

À 

| 

| 

| 

d 

| 


Fig. & — Resposta à função escalão. Gráfico obtido no registador 

X Y segundo o esquema da Fig. 7. Encontram-se assinalados os 

valores utilizados no cálculo para o caso indicado na 3.º linha 
do Quadro | 


A partir dos valores experimentais dos elementos 
do filtro calcularam-se os polos e zeros de Hlis]/s e a 
partir de colecções de valores de u; (t) quer obtidas com 
o voltímetro quer obtidas do gráfico do registador uti- 
lizou-se o programa descrito para obtenção dos mesmos 
valores, 
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No Quadro | encontram-se alguns dos resultados 
obtidos, apresentando-se na figura 9 os resultados im- 
pressos pelo programa para um dos casos considerados 
(2.º linha do Quadro |). 


QUADRO 1 


PROCESSO 
UTILIZA DO 


Cáleulo a parir 
dos valores ex- 
perimentais dos 
clementos 


Frograma-Volt, 


gm cAt=60 pr - 
|- 16 pontos 
- Grau = 8 


Frograma-CGrd- 


Os resultados apresentados na figura 9 referem-se 
a um caso em que a análise foi feita com ordem supe- 
rior em duas unidades à do sistema em estudo, Tal como 
foi citado na descrição do programa, aparecem duas das 
constantes (sublinhados a traço interrompido) cujos mó- 
dulos são bastantes menores que os módulos das restantes 
(cerca de 5% do menor módulo das outras constantes), 

Como pode verificar-se o polinómio numerador a 
que se chega, não é real, sendo as partes imaginárias dos 
coeficientes devidas aos erros cometidos na execução 
dos cálculos (valores da ordem de 10-4 da parte real). 

Neste exemplo, os zeros do númerador deveriam 
astar todos no infinito, Tal não sucede no cálculo devido 
aos erros referidos e ainda aos erros na medida do 
tempo em relação à aplicação do escalão. No entanto, 
observa-sa que todos eles têm um módulo maior que o 
módulo de qualquer polo (Fig. 10). O pequeno valor 
dos coeficientes dos polos resultantes do aumento da 
ordem da análise traduz-se pelo aparecimento de zeros 
junto desses polos, 

Em relação aos valores do Quadro | verifica-se que 
os polos obtidos pelo programa concordam com as 
obtidos por cálculo a partir dos valores experimentais 
dos parâmetros, a = 2% do maior módulo. 


4. CONCLUSÕES 


Mostrou-se que a partir dos valores da resposta 
impulsiva tomados em, pelo menos, 2 g, instantes 
separados de At é possivel identificar completamente a 
função caracteristica dum sistema linear de ordem 9p- 

Para isso estabeleceu-se a equação de recorrência 
a que esses valores devem satisfazer e relacionou-se 
a solução geral desta equação com a expressão geral 
da resposta impulsiva dum sistema linear. 

Elaborou-se um programa para o cálculo automá- 
tico dos polos e zeros da função característica a partir 
dos valores das observações da resposta impulsiva, O 
qual foi ensaiado com exemplos numéricos e experi- 
mentais, tendo-se obtido resultados satisfatórios, 
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dicado na 2,º linha do Quadro | 


" 
: PIBUAT =» 
ed + o F 
: ' Poa ta 
ã Pi "fe dao th 
é Pica 4 + 
E 
F 
ps R E o a L 
+ UN EA 1% 
E NÉ E T a To) TT e À f i “4 f E jr ET. 1 SEE Er ai 
POA VU SUA LO ILPELELCO d TOU Ir iys4sTO-" DU ERICO) 1994" 0=) 
a MRE, MC re SD Ta CO CR 
e ; À E q ak 7 o nr ê É FR is vi r 7 Fo : i o é 
CLS LU AUS cteve gy) 4 A É do é a ICO “00 JTIQSES*U=) SO JELG9RAS*D= "20 IT92BDT' j 
Bis Lt atira 
| E; b À LM 4 ; l = 1 z F ) r a Ud : z P 
(ALUS ICE Tra E] L Rm e dom o | 1 | [Ay [URGE E P fla] [9L= dE EH Ru , t ga DR A a [64º ud | USO = Sp o) E de = * k "= 
[1 - dll NS ip= E) RR" S | a 1 mid | rd i=+t no ] | Eli=3+ ] a E ' hd E fa É | f e RR, 4 é 
4 Ly O F, L L Á ' j 4 t+ | É À Ê te PA te ) NI Is ELE li | l Rue so Es | ta (1 E Í 
É se Li Lg ; =] e err a Ar | ; . ] i 
|? ” ?1 E As ais “ ea Re o EE prod po | TE TIN Ml | 
a 
e E A E E SE SS O SR IS CLESNS a EO IS OLIGELS Et CABLIAGis to 2 DONE E ESCILCG ES KO E io ee EA E Ss E SS pis gia e rir Sie ST E iii pg pa e ç 
LO ] | = 1 3 | E 1 Fi - 7 1 ja k 
ulo Clu-9LeR U= + 4 Er 14590 ) 400 4 *u=* gs: “U- H10€902*0 
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Não se considerou o problema da escolha do in- 
tervalo 4t e do número de observações a tomar para 
obter a melhor precisão dos resultados para uma dada 
imprecisão das medidas, 


Fig. 10 — Localização dos polos e zeros obtidos pelo programa 
para o caso indicado na 2.º linha do Quadro | 
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CONVERSOR DIGITAL ANALOGO TEMPORAL 


RESUMO 


Ensaia-se a sintese de um pequeno computador 
especialmente projectado para controlar os movimentos 
de uma antena helicoidal, destinada à captação de foto- 
grafias meteorológicas enviadas pelo satélite ESSA - 8». 


NOTA PRÉVIA — Para comando automático das 
antenas usadas na recepção de sinais de satélites arti- 
ficiais, no CENTRO ESPACIAL DA MULEMBA, fazia 
parte dos programas para 1973, o projecto de uma 
unidade digital, cuja construção seria levada a cabo 
nos próprios laboratórios do Centro. 


Acreditando, todavia, que esta tarefa poderia ser- 
vir de boa prática e estimulo aos alunos do Grupo de 
Estudos de Cibernética da Universidade de Luanda, foi 
o director do C. E. M. contactado afim de lhes forne- 
cer os elementos-base para o referido projecto. 


O presente trabalho, executado por Rodrigo Mar- 
tins e Isaura Simões, representa assim, simultâneamente, 
uma valiosa colaboração para o Centro, e excelente 
prática de estudo para aqueles alunos da Universidade 
de Luanda, 


1 — INTRODUÇÃO 
1.1 — Objectivo e finalidade do programador 


O programador que nos propomos construir, serve 
para comandar automáticamente uma antena do tipo 
helicoidal. Os dados (altura e azimute) são fornecidos 
ao programador, minuto a minuto, durante o tempo má- 
ximo de 25 minutos. Uma vez que são fornecidos com 
bastante tempo de antecedência, podemos arranjar um 
sistema, que armazene o fluxo de informação, de modo 
a ser removivel sempre que se deseje, 


Assim, podemos apresentar desde já, as seguintes 
questões: 


a) — Como são fornecidos os dados? 


b) — Qual o processo a utilizar para introdução dos 
dados na memória? 
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RODRIGO MARTINS 
ISAURA MARTINS 
Grupo de Estudos de Clbernética 


Universidade de Luanda 


SYNOPSIS 


alt is diseussed the design of a special purpose 
computer, capable to control two movements of an 
helycoidal antena which receives meteorological photo- 
graphies from ESSA -B satelyte.» 


c) — Qual a «lógica» a usar para o armazenamento 
dos dados. 


d) — Como se há-de «tirar» a informação da me- 
mória. 


e) — De que modo a informação «atirada» vai actuar 
sobre o «sistema de saida». 


As respostas, a estas questões, definirão o modo 
de proceder para a construção do programador. 


a) — Os dados são fornecidos na forma décimal, 
sendo a altura representada por dois digitos 
e o azimute por três digitos. A introdução 
destes dados, far-se-á por meio de uma uni- 
dade de entrada a designar. 


b) — À introdução dos dados na memária será feita 
por um selector de canais; Antes, deve-se 
fazer um armazenamento prévio destes, de 
modo a que não haja perda de informação. 
Este armazenamento prévio será feito no 
«BUFFER» de entrada. 


ci — O tipo de memória que vamos utilizar é cons- 
tituido por núcleos de ferrite os quais arma- 
zenam a informação sob a forma octal, pelo 
que, dever-se-á construir um conversor déci- 
mal-binário /octal-binário antes do «BUFFER» 
de modo a que as saídas seleccionadas já se 
encontrem na forma octal-binária. 


dj) — O fluxo de informação sai da memória por 
meio de um selector de canais. 


el -— A informação, uma vez tirada da memária 
deve ser também armazenada num «BUFFER» 
para que haja uma maior estabilidada de 
circuito. 
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f|l — À informação, uma vez nos registadores de 
saida, vai actuar sobre um comando de «em- 
braiagem» da antena helicoidal, durante uma 
quantidade de tempo, dependente do valor 
lógico da informação transmitida. 


Nestas condições estamos aptos a apresentar a 
arquitectura do programador que nos propomos cons- 
truir, 


1.2 — Resumo da sua arquitectura 


As unidades principais, constituintes do programa- 
dor, foram já focadas no parágrafo anterior, Essas uni- 
dades, são indicadas por meio do esquema de blocos 
da fig. 1-a: 


— UNIDADE DE ENTRADA — Serve para a recep- 
ção dos dados e é constituída por «THUMBWELL SWI- 
TCHES+». Incluída na unidade de entrada temos um sis- 
tema «apontador» dos dados iniciais (para um primeiro 
controle) constituido por contadores, «Drivers» e «Ni- 
xies». 


— CONVERSOR  BINÁRIO  DECIMAL - BINÁRIO 
OCTAL — Esta unidade serve para fazer a descodifica- 
ção da entrada binária-decimal na lógica em qua é 
actuada a memória. 


«BUFFER» DE ENTRADA — Sistema que serve 
para se fazer o armazenamento prévio dos dados em 
octal-binário, antes destes serem colocados nos núcleos 
da memória. 


— MATRIZ DISTRIBUIDORA (OU COMUTADORA) 
DOS DADOS DE ENTRADA — Sistema utilizado para 
a colocação dos dados de entrada nas respectivas célu- 
las de memória (matriz de diodos), 


- MEMÓRIA — Armazém de dados constituido 
por núcleos de ferrite. 


— MATRIZ SELECTORA DE CANAIS — Sistema 
utilizado para «retirary a informação da memória e colo- 
cá-la no «buffer» de saida (matriz de diodos). 


- «BUFFER» DE SAÍDA — Unidade que serve para 
reter temporariamente a informação, antes desta passar 
para a unidade de saída. 


— UNIDADE DE CONTROLE — —Esta unidade 
astã colocada a seguir ao «Buffer» de saida, e é 
constituida por um sistema tipo agulheiro, que permite 
a passagem da informação para a unidade seguinte, 
minuto a minuto. 


- CONVERSOR DIGITAL - ANALÓGICO - TEMPO- 
RAL — Sistema, que vai converter a informação lógica 
de saida em: 
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1.º — diversos niveis de tensão. 
2.º — diversos intervalos de tempo (os níveis deter- 
minados anteriormente). 


Anexo a esta unidade, encontra-se um sistema adi- 
cionador temporal (monostável e porta OU). Todo este 
conjunto, vai actuar sobre a embraiagem da comando 
do motor de antena. 


2 — ESTRUTURA DO PROGRAMADOR 
2.1 — Opções tomadas, sua explicação 


a) — Para UNIDADE DE ENTRADA, escolheu-se 
um «Thumbwell Switcher, pois este, dá-nos 
directamente a informação décimal, na forma 
binária-décimal. 

Também se poderia ter escolhido para 
unidade de entrada um sistema contador, 
accionado por relés (tipo máquina de escre- 
ver), com perfeita adaptação ao sistema es- 
colhido. Contudo, optamos pelo «thumbwell 
switche» pelo seu reduzido tamanho, boas 
características ce funcionamento, custo rela- 
tivamente baixo e alta precisão. 

Contudo este sistema apresenta um in- 
conveniente, no seu manuseamento, pois o 
modo como são inseridos os digitos decimais 
não é nada cómodo, 

Este defeito, não existe no sistema de 
relés. Contudo, o sistema de relés não apre- 
senta uma acuidade tão grande, devido à 
inexistênsia de um memorizador temporário, 
pelo que, a solução ideal seria aplicar o sis- 
tema interno do «thumbwell switches» a um 
comando tipo «máquina de escrever», O qua, 
posteriormente pode ser feito. 

Esta unidade encontra-se representada 
na fig. 2. 


bj — Para conversor-binário-decimal/binário-octal 
utilizamos um sistema, tal como se mostra na 
fig. 3 ja, b), o qual é bastante complexo. 

Juntamente com este sistema, apresen- 
tamos um outro, o qual apresenta grandes 
simplificações em relação ao sistema ante- 
rior. 

A entrada para este sistema é feita em 
séries, uma vez que a saída do T. S. é feita 
dessa maneira, À saida desta unidade é que 
pode ser em paralelo ou série, consoante o 
tipo de conversor utilizado (3a ou 3b), 


c) — O «BUFFER» DE ENTRADA, é constituído 
por um registador de 12 bits, cuja entrada é 
feita em séria ou paralelo, consoante o tipo 
de conversor escolhido, Contudo, optamos 
pelo segundo sistema do conversor, por ser 
bastante simples, pelo que, podemos dizer 
que a entrada para os registadores é feita 
em série. 
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A saída dos registadores é feita em 
paralelo pois interessa-nos tal tipo de saida, 
para que a unidade comutadora ou distribui- 
dora de canais possa abranger todas as cé- 
lulas da memória escolhida. 


dj) — A MATRIZ DISTRIBUIDORA (ou comutadora) 


é constituída por diodos e possui 12 entra- 
das e 4096 saídas. Optou-se por este tipo 
de matriz por ser relativamente barata com- 
parada por exemplo, com um contador de 
4000 palavras. Assim, para este caso, a co- 
mutação é feita directamente, no caso de 
se utilizar um contador de 4000 bits (adress- 
-register) teriamos de fazer um programa 
prévio para a introdução dos dados, na me- 
mória. Na solução escolhida, os dados são 
introduzidos pela mesma ordem de re- 
cepção. Ássim, dever-se-á fazer uma «tábua» 
onde figurem as «posições» de ligação da 
matriz de diodos (saídas) aos respectivos 
núcleos de memória, de modo a saber-sa 
onde serão armazenados os «bits» duma in- 
formação,, para se conhecer qual deve ser a 
sequência de saida da informação. 

Tal tábua será feita mais tarde, depois 
do teste de todos os constituintes. 


e) — MEMÓRIA — Utilizou-se uma memória de nú- 


cleos de ferrite, com capacidade de armaze- 
namento de 4000 bits, com as seguintes ca- 
racterísticas: 


Corrente de magnetização dos núcleos 500- 
-600 mA; 

Tensão de entrada — 30-60 mvV; 

Núcleos dispostos numa matriz de 80 x BO. 


f)— MATRIZ SELECTORA DE CANAIS — Utili- 


zamos a mesma matriz de entrada, com um 
funcionamento inverso, de modo a que a in- 
formação saia em série. Exista à saída um 
«DELAY» que só deixa passar um bit de cada 
vez (nas 12 saídas), de modo que a infor- 
mação sala em série, e não em paralelo. 

Também poderiamos tirar a informação 
em paralelo, originando assim uma modifi- 
cação no registador de saída. 


9) — «BUFFER DE SAIDA» — É contituido por re- 


gistadores (F. F.), os quais possuem uma 
única entrada e 24 saídas, (2 registadores 
de 12 bits em série) no caso de se optar 
pelo primeiro tipo de matriz selectora; ou 
possui 12 entradas em 24 saidas (2 regis- 
tadores de 12 bits associados em paralelo) 
no caso de se optar pela segunda solução 
da matriz selectora. 


|) — RELÓGIO CONTROLADOR — o funcionamen- 


to da porta de controle, quer dos F. F. selec- 
tores, e dos registadores é feito por um re- 
lógio que produz um nível «alto» de minuto 
a minuto (tempo de espaçamento entre duas 
sequências imediatas) com uma largura À 1. 
(20 us por exemplo). 

Assim, este impulso é aplicado directa- 
mente ao C., P. dos registadores, podendo-se 
desligar o da entrada ou da saida consoante 
o fluxo, está a entrar na memória ou a sair 
desta. O impulso de saída é accionado pelo 
comando dos canais. Assim, quando pre- 
tendemos retirar o fluxo de informação 
que se encontra armazenado na memória, car- 
regamos no interruptor dum dos canais, intro- 
duzindo o impulso de relógio que vai actuar 
nos F, F. de minuto a minuto, durante um 
tempo At; por sua vez, existe um monostável 
com um tempo 12 At, que vai fazer o RE- 
SET do interruptor provocando o acciona- 
mento do S. e assim sucessivamente, 

O contador que vai permitir a abertura 
das portas do relógio tem um tempo de co- 
mutação de 24t, Este contador, pode ser 
substituido por um monoestável com um 
tempo de actuação de 2At, de minuto a mi- 
nuto. 

Só nos interessa que o relógio provoque o 
accionamento do sistema programador du- 
rante 25 minutos de accionamento, o que 
corresponde à leitura de 600 bits (altura 
de azimute) correspondentes a 2 canais, 
e mais de 12 bits do canal seguinte, 

No fim desta tempo (25A4t) o relógio pára, 


impedindo a perca de qualquer informação que 


se encontre armazenada nos núcleos seguin- 
tes, 


ii — FONTE DE ALIMENTAÇÃO — A polarização 


da matriz de diodos e do conversor A/T, é 
feito por uma pilha de +15V. A partir desta 
fonta por meio de resistência divisoras de 
tensão, pode fazer-se com que a tensão de 
cada núcleo seja de facto mantida a 30 mv, 
Este aspecto da fonte vai ser tomada em 
conta durante os ensaios experimentais, 
pois agora, estamos só interessados em 
determinar a «estrutura e arquitecturas do 
programador prôpriamente dito, tal como se 
vê na figura nr, 1, 

Esta fonte de alimentação é a mesma para 
oT.s. 


2.2 — Modo de proceder para armanezar ou retirar 


o fluxo de informação. 


Suponhamos que foi introduzido o número 56 na 

hj) — À UNIDADE DE CONTROLE escolhida é do unidade de entrada, sendo o algarismo 5 introduzido na 

tipo «agulheira» e de bastante fácil manusea- parte correspondente às dezenas e algarismo 6 às uni- 
mento. dades. 
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A entrada de cada um dos contadores aprece-nos 
a seguinte expressão: 

0101 (nas dezenas); 0110 (nas unidades). 

A informação sob esta forma, passa para o con- 
varsor decimal-binário, octal-binário, aparecendo-nos 
à saida do conversor a seguinte sucessão. 


0000 0111 0000 
= —— 
MSB LSB 


Nota: MSB — Bit mais significativo; 
LSB — Bit menos significativo. 


Esta sucessão é armanezada no registador do LSB 
a MSB. 

Consideremos agora que a saida S, da matriz de 
diodos está ligada à unidade de memória correspon- 
dente a X;; Yi; S:a X:,Yi;; S/a X,Yi:;S aX, Yi Sa 
Ki Yi; Se a Xe Ti Sa Me Vijcccaccs Su O Xis Vis 

Nestas condições, quando houver mais um impul- 
so no CP. do registador ou quando aparecer um novo 
fluxo de informação, a informação anterior será armaze- 
nada nos núcleos correspondentes à 1.º coluna da me- 
mória e às doze primeiras linhas, ficando armazenado na 
memória os seguintes bits: 


Para se retirar a informação 
Y, da memória, carrega-se no «inter- 
ruptor» do canal 1, que abre uma 
porta e permite a passagem do 
impulso de relógio que vai actuar 
sobre um F.F.D. (CP-=-D) com 
um tempo de actuação At. Ássim, 
passado um tempo dado (próximo 
de minuto), o flip-flop actua sobre 
os 12 primeiros núcleos excitados 
pela saída 5, e, como matriz de 
diodos agora está ppolarizada in- 
versamente, a informação flui atra- 
vês desta, dando uma saída em 
séria ou em paralelo, que vai para 
um registador de saída. 


M.s.B. 


00004 +4+0 0000 


L.5.B. 


Após o instante At o impulso em D 

vai a zero, mantendo-se neste estado duran- 

te At, fazendo com que 5 — 1, e, como 3 

está ligado a Y, e não existe aí nada arma- 

zenado vai-nos dar uma saída nula e assim sucessiva- 

mente. 

Nestas condições, no registador de saida, e nos últi- 

mos 12, ou nos 12 primeiros FF-3 aparecer-nos-ia a 
informação desejada. 


0000 0000 0000 0000 0111 0000 
ou 

0000 0111 0000 0000 0000 0000 

MSB LSB 


Como passado um minuto, um monostável (ou 
contador de 12 bits) actua sobre as portas existentes 
à saida do registador, esta informação passa para O 
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conversador D/A de modo a que o bits M.S.B, vá corras- 
ponder a um nivel baixo i. e., as entradas ficam inver- 


tidas. 
0900 1110 0000 0000 0000 0000 
Vai para o conversor 1 Vai para o conversor 2 
(Altura) (Azimute) 


À saída do conversor 1 vamos ter uma tensão que 
nos é dada por: (Ver o correspondente ao circuito 


D/A-T) 
Voa — Vint 
1 1 1 1 
[1 4 ! 4 8 s 4 
1 1 1 1 
fe Bj do e! aid ud 
1 1 o 
ea D, + ET Dio + 3 D; + “gs Du) Va 


D,=D,.= D,=D,=D.=D=Daz= Dun=DHz=zO0 


1 1 1 
“ig A o o E) a 
VR=10 Vou 5V 
0,875 
os MA 


Voa | =27,343 x 10-* V 
Voa, = 13,6715 x 10—" V 
3 


T; = 0,73 s 
Tá 1,46 s CM, 


lt 


Se a tensão à saida do registador fôr de 10 V, a 
informação a actuar durante 0,73 s + o tempo de acrés- 
cimo mínimo dado por um monostável sobre a embraia- 
gem do motor ocasionando o deslocamento do mesmo, 
Se Va =5 V, o tempo de actuação sobe para 1,46 
s + Tmonos x 

Assim, no fim deste ciclo, houve a tansformação 
da «informação digitos, em «informação tempo» que 
actuou sobre o motor provocando o deslocamento em 
altura deste. 

Como não foi introduzido qualquer azimute 
Va — O não existe tesão à entrada do conversor dois, 
pelo que não houve deslocamento ocasionado por esta 
variável de entrada, 

Todos os outros ciclos do fluxo de informação, 
processar-se-iam de um modo análogo ao indicado. 
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2.3 — DESCRIÇÃO DE CADA UMA DAS UNIDADES 
QUE CONSTITUEM O PROGRAMADOR 


2.3.1 — UNIDADE DE ENTRADA 


Como já se disse anteriormente, esta unidade é 
constituída por um «thumbwel| Switchs 1248 N com 
uma entrada para as unidades de centenas; unidadas 
de dezena e unidade de unidade (ver fig. 1). 


A tábua de verdade é a seguinte: 


h3 
a | 
Co 


INDEX f 
2 0 o 0 o 
1 || q o 0 0 
2 0 1 oo 
3 1 1 o q 
4 0 o 1 0 
5 1 0 1 q) 
6 o) 1 1 1) 
7 1 1 1 0 
8 o o o 1 
9 1 o 0 1 


Para que haja um prévio controla dos dados que 
estão a ser inseridos no programador existe uma unidade 
ufixadora temporáriandos bits à entrada, contituida por 
3 contadores (um para cada unidade) SN 7490; 3 
adrivers» SN 74414 (um para cada unidade) e 3 enixies» 
(uma para cada unidade). 

A informação, sai em séria do T.5S. e, sobre a forma 
binário-décimal; entra sob esta forma para o conversor 
binário-décimal/binário-octal. 

Por que se proferiu esta unidade já foi explicado 
anteriormente, e, como se trata de uma unidade já 
«contituída» nada mais acrescentaremos a seu respeito. 


2.3.2 — CONVERSOR BINARIO-DECIMAL / BINÁ- 
RIO-OCTAL 


À existência deste conversor, deve-se ao facto 
da memória que vamos utilizar armanezar os bits sob 
a forma octal. 


MULTIPULICADOR/TO 


cida 4 16 6 4 2 1 
1/0001 0001010 
2 106040 o 0 o 100 
3/1/0011 o o 1147) 
410100 0101006 
5/0101 0119066016 
610110 114141 0:00 
7) 61144 1000110 
8/1200 11004060 % 
& | 7001 Ot Aa 
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512 256 128 64 32 16 8 4 


Como já se disse, anteriormente, este conversor 
pode ser construido de dois modos diferentes. À primei- 
ra solução é a apresentada na fig. 2-a, em que se indica 
o esquema de blocos como também o número de uni- 
dades necessárias para a sua constituição, 

À outra solução baseia-se na anterior com a van- 
tagem de nos trazer simplificação de blocos (fig. 2-b). 
Baseando-se no esquema de blocos anterior, vamos des- 
crever o funcionamento do conversor. 

A entrada do conversor é feita por contadores de 
décadas (3 correspondendo às 3 entradas — U. — D, — 
— C.) os quais possuem um circuito inibidor para eli- 
minar qualquer ocorrência de algarismos incorrecta. Es- 
tes também servem de «armazenadores» prévios. 

O circuito conversor em si, é constituido por cir- 
cuitos somadores e multiplicadores por 10 e 100; sendo 
estes últmos constituidos por «Shift-left Register'sa. 

Na folha em que se indica o esquema de blocos 
inicial, encontra-se uma tábua de verdade para as uni- 
dades multiplicadoras, 

As três saidas do T. S. (V.D.C.) são bifurcadas à 
entrada do conversor, entrando cada saida para um con- 
tador de décadas e para um «SHIFT-LEFT Register de 
2 bits»; excepto a combinação de bits correspondentes 
às unidades, que só entram para um contador de dé- 
cadas. 

Quer a entrada do contador quer a do registador 
é feita em série D, a C.). 

Assim, a combinação de bits correspondentes às 
centenas entra para o contador o qual os memoriza e 
os reproduz à saida passado um instante At. (4 saídas 
em paralelo); por outro lado, o registador ao receber 
a informação, (4 bits) por meio de dois impulsos à entra- 
da (CP) provoca o deslocamento para a esquerda da in- 
formação, ocasionando o aparecimento da mais dois 
zeros à direita do número (registador de 6 bits). 

Pelo mesmo comando (CP) a informação deixa o 
registador por uma única saída, (série) a qual sa vai ligar 
à entrada de um F.A. 

Como a saída do contador á um paralelo, neces- 
sitamos de um conversor paralelo-séris de modo a que 
a entrada da informação no F.A. seja feita em série. 

Este conversor (ver fig. 3) é constituido por um 
registador de deslocamento cuja função de controls é 
realizada por um contador em anel. 


MULTIPULICADOR/ 100 


o OQ 1 o 0 O 

[ê, 1 1 o O 1 o 0 OQ 
1 o O 1 O 1 O O 
1 1 o 0 1 oO 0 0OoãÃ)lf 
1 1 1 1 1 O 1 o OQ 
o O 1 O 1 1 O 0 OQ 
o 1 O 1 1 1 1 o O 
1 o O 1 o o À0Oo go 
1 1 o 0 oo 1 o O 
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O registador de deslocamento, que compreende 
uma entrada série Is (directamente ligada à entrada J 
ou D do primeiro F.F.) mormalmente conservada no 
valor lógico O, é constituida por 6 elementos, sendo os 
primeiros 4 destinados ao armazenamento da informação 
paralelo, o quinto destinado a controle e o último a uma 
linha de atraso de 1 tempo. 

O contador de anel, também constituído por 6 ele- 
mentos (F.F.) é tipo «self-starting» and self correcting 
ring, «counter», — o qual dispõe de uma entrada de 
«direct clear» Ry comandada pelo conjuto de controle. 
O gerador de impulsos é constituído por um oscilador 
de cristal de quartzo seguido de um «Schmitt Trigger» 
(dispara só no nível alto) e fornece uma sequência con- 
tinua de impulsos retangulares a intervalos de tempos 
iguais, cada um dos quais constitui uma unidade de 
de tempo, ou seja | tempo característico do sistema. 

O conjunto de controle, compreende um F.F. R.5. 
básico, dois circuitos lógicos «AND» designados pelos 
nº 1 e 2, um inversor | de entrada S, e de saída di- 
rectamente ligada à entrada S! do quinto F.F. do regis- 
tador de deslocamento, 

Os circuitos «AND», designados por D, €, B, À, 
têm as suas saidas directamente ligadas às entradas S 
dos primeiros 4 F.F. do registador. 

Para a análise do funcionamento do conversor pa- 
ralelo-séria acima descrito partiremos dum estado inicial 
caracterizado pelas condições 0 = O, S, Ip =le=— 1, 

-| — O e pela «limpeza» total do registador (So = 
=D=C= B=-A=O). Porser R, — O, a condi- 
ção O — O acarreta R, = O o que «limpa» totalmente 
o contador em anel, que apresentará portanto o valor 
lógico O em todos os seus 6 F.F., como se indica na 
fig. À mesma condição O = O bloqueia a transferência 
para C.P,, através do circuito lógico «AND» número 2, 


dos impulsos gerados por D,. À condição S, — O 
implica S, = | na saída do inversor |; no entanto esta 
condição S, = | não é trasferida, através do circuito 


lógico«AND» número | para a entrada S do F.F, RS, de- 
vido ao bloqueio exercido pela condição So = O, o que 
implica 5 = O neste F.F. Finalmente, como em R.5, 
sa tem R no mesmo estado lógico que no último F.F. 
do contador em anel, tem-se ainda R = O no estado 
inicial do conversor, 

Averiguemos como se comporta o conversor para- 
lelo-séris quando na sua entrada paralelo se aplica uma 
informação contituida pela série de bits 0101 equivalente 
ao algarismo 5. 

Admitamos que a entrada desta informação paralelo 
se realiza premindo um botão não representado na 
fig. 3, |. e que existe informação paralelo à saída do 
nosso contador por se ter aplicado um 5, na «unidade 
de entrada» correspondente às centenas. 


Tal informação, aos terminais do contador faz com 
que S, =lelp la ;ly ;l, = 0101, condições que 
se matêm nas entradas, apenas enquanto o contador 
não comutar de estado e regressam a zero logo que 
esta se dê, Então, a presença da informação aos termi- 
nais do contador faz So = 1 e DCBA = 0101 como se 
indica na figura. Enquanto existir o fluxo da informação, 
tem-se S, =1>5, =0O e portanto a condição So 
= 1 não é, por enquanto, transferido através da porta 
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«AND» n.º 1, para entrada S de RS, mantendo-se S 
=0e 0-0. 

Ao dar-se a comutação no contador, os valores lógi- 
cos de S,, lj, 1,, Ig. |, regressam a zero mas manten- 
do-se So = 1 e DCBA = 0101 no registador. 

Agora as condiçãos So = 1 de entrada da porta 
«AND» n.º 1 fazem com que S = 1 em F.F. R.5. 

Esta condição acarreta a mudança de estado de 
dae RS de 0 = 0 => 0=71,0 que acarreta o desbloquea- 
mento da porta «AND» n.º 2. Nestas condições, come- 
çam a ser aplicados impulsos CP ao registador e ao 
contador em anel. 

A aplicação destes impulsos em CP faz deslocar 
para a esquerda, instante a instante, no registador, a in- 
formação que vai ser transmitida, pela saida, na forma 
de informação série lo, precedida ou não do bits 1 de 
controle representado à tracejado, que o nosso caso é 
inibido à saida. 

Às eventuais mudanças do valor lógico So resul- 
tantes deste deslocamento da informação acarretam 
mudanças identicas em S, que não fazem mudar de 
estado RS enquanto = O. Por outro lado devido ao 
facto de sar IS — O os impulsos aplicados ao registador 
vão, passo a passo (tempo a tempo, substituindo 
por «Zeros» os valores lógicos inicialmente regis- 
tados. Nestas condições, depois de terminado o 6.º im- 
pulso, toda a informação foi transmitida pela saida do 
registador na forma série e este apresenta-se limpo. 
À aplicação do primeiro impulso no contador em anel 
faz aparecer o valor lógico 1 no 1.º F.F., valor esse que 
se vai deslocando instante a instante da direita para 
a esquerda. 

Nestas condições, ao terminar o 6.º impulso,o ulti- 
mo F.F. do contador apresenta o valor lógico 1. Este 
facto faz com que em R.5. onde nesta altura já se tem 
Ss - O, seja R — 1. Nestas condições dá-se, em RS, o 
retorno de O ao seu valor inicial (0=1>0=0). 

Este retorno faz bloquear o circuito «AND» n.º 2, 
O que acarreta a cessação do fornecimento de impulsos 
a CP e origina R, — Q = O, condição que limpa o 
contador em anel. Esta análise mostra que, depois da 
da transmissão da informação da forma sérias, o convear- 
sor fica automáticamente preparado para nova conversão 
paralelo-série. 

Assim a informação é transmitida por meio de duas 
entradas (A; B,) ao circuito somador. Uma delas 
está directamente ligada ao registador de deslocamento 
de dois bits e a outra ao conversor paralelo-séria, 
Como queremos que a cadência de entradas no circuito 
somador seja a mesma das duas linhas, poderemos uti- 
lizar o mesmo oscilador (e ring counter) quer para o 
conversor quer para o registador de deslocameneo de 
dois bits. O F. A. possui 3 entradas e duas saidas 
estando duas das entradas já ocupadas e a outra vai 
ligar à entrada de um circuito delay (C,), com o outro 
terminal ligado à saída 5, do F. A. para onde vai o bit 
de transporte, sendo a outra saída S, a saída própria- 
mente dita do F. A. 


O circuito delay apresenta o atraso de 1 bit. Assim 
se em À, tivermos a sucessão 0101 e em B, 110100 


ns 


deslocamentos de dois bits. 
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ão 1.º impulso em A, entra 1 e em B, O saindo 
em S,=1 e S;= 0. Como o circuito de delay está 
regulado para ter o atraso de um bit, o valor da saída S; 
vai-se somar à nossa sequência de entrada. 

Para o 2.º impulso tem-se A, = 0; B,=0;C,=090; 
dando em S,=0 e 5,=0. 

Para o 3.º impulso tem-se A, = 1; B,=1;C,=0. 


Agora, na saida S, temos: S,= 0 e S;= 1. Como 
o circuito de delay introduz um atraso de um bit, para 
a próxima sequência já se tem C; = 1, Assim, para o 
4.º impulso temse A,=0;B,=0;C=1=—> 5S=1, 
S,= 0. Para o 5.º e 6.º impulsos a saída S,= B, e 
S, = O constante por já não existir informação em As. 

Nestas condições, à saída deste somador teriamos 
a sequência 0171001. 

Esta saída liga a um registador de deslocamento 
de um bit-séria. À saída deste registador teremos a se- 
quência 0110010, a qual vai entrar para outra unidade 
de F, À. com delay idêntica à anterior e para um regis- 
tador de deslocamento para a esquerda de 2-bits, apa- 
recendo à saída deste último a sequência 011001000. 
A saida deste registador constitui a nova entrada para 
o F.A. anterior. Por um sistema idêntico ao anterior, 
teriamos a saída deste somador a sequência 011111010 
que iria entrar num registador de deslocamento de um 
bit, aparecendo à saida deste a sucessão 0111110100 
a que vai ligar ao somador 5, como se vê na figura 
n.º 2b. 

Se suposermos que o algarismo das dezenas é 9, 
teremos também à entrada, um contador e um circuito 
conversor paralelo-séria acopolado a este e um regista- 
dor de 2-bits cujas saídas vão ligar a um circuito soma- 
dor com DELAY, que nos vai dar uma saída S, que 
apresentará a seguinte sequência 101101. Esta saída 
vai por sua vez ligar a um registador também séris de 
deslocamento de um bit, cuja saída se liga a um circuito 
somador 4 como se vê na Fig. 2b. A este circuito 
somador vai ligar-se a saida do conversor paralelo-sério 
do contador das unidades. 


Suponhemos que o algarismo das unidades é 8. 
Nestas condições à entrada do circuito somador teria- 
mos as sequências 1011010 e 1000 que iriam originar 
uma saída 1100010. 

A saída deste somador liga à entrada do somador 
5 onde val ter também a sequência, já determinada das 
centenas, 0111110100. 

A esta saida podemos ligar um nivel fictício sem- 
pre ao nivel zero, pelo que, podemos dizer que o alga- 
rismo apresenta-se à saída com a forma 


001 001 010 110 


massas 


Na forma octal temos o número (1126), — (598)... 
Se compararmos este circuito com o anterior, Fig. 2a, 
verificamos que houva uma grande melhoria pelo que, 
resolvemos chamar a este sistema de semi-hardware/ 
/semi-software., 
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2.3.3 — Registador de entrada 


Esta unidade, tem como função, armazenar os da- 
dos provenientes do conversor binário-décimal/binário- 
-potal, A informação proveniente do conversor é sem- 
pre constituída por 12 bits (ou 9 bits com + 3 níveis 
zeros de modo a uniformizar a sequência) série que 
antram no nosso registador de deslocamento com saldas 
am paralelo, 

Este registador (Fig. 6) possui um relógio qua co- 
manda as sequências de saída e entrada de dados de 
modo a que as sequências sejam sincronas. 

O registador, é constituido por 13 F. F. que cons- 
tituam o circuito integrado. 

O motivo que nos levou a optar por um sistema que 
nos desse as saídas em paralelo já foi considerado, e, 
o processo pelo qual se faz o armazenamento é análogo 
ao indicado para o conversor paralelo-séria, agora des- 
crito em sentido contrário. 

Assim, se todos os F.F, possuírem um dada in- 
formação, a nova sequência faz com que esta transite 
para a matriz distribuidora, ocasionando um «vazio» no 
registador que vai ser preenchido com a nova informa- 
ção. Este processo repete-se ciclicamente. Todo este 
controla é feito pelo CP de registador. Para o nosso 
caso teriamos: 


CP Data A Ê D Ê a + J k 

- - Ú Ú l Ú Ú Q q ' A [a L Ú L 
k l 1º [E O Ú Ú Q Q ú a O U E E 
2 bl bl 1a O oe Ooo e o c 
3 bz ba 1 1 q Ú O O Ú o Fr Q LP) o 
4 bj t3 &b LL 1º 0 O Ooo o og o 
E) ba b4 b t2 bl 1º O [8 Õ O O [à O 
b b5 t bs bl b? b 1* O Ú t 1 ( Ú Ú 
7 bã k . b b bL 1º 4 [à é E 
8 b7 b7 bt b) bã bl) t bi a 6 EU [a O ; 
E] bã aa b7 bt > bã t b bl 1% O Ú L 3 
10 b9 bg A b7 tê b5 tb b3 bz bl “ 
Lil b1LO bH by ba b7 tb b bt bã E& b ] O 
Le bill tbil tblo bo bB b7 bh bh b4 bã b?2 bl WU 
13 b12 bl? bill blO b5G b& b7 tê bS MM bi ba bi 


Em que 1* é o impulso de comando que vai «abrir» as 
saídas (ver Fig. 6) e N + é o F.F. de controle. Para 
os registadores de deslocamento de 1 ou 2 bits estu- 
dados anteriormente, o bit de controle era zero com 
saída própria e um inversor na sequência de saídas i. e. 
de modo a ter-se a sequência inicial, mas deslocada 
para a esquerda com uma única saída. 


2.3.4 — Registador de saida 


Como a nossa matriz selactora (ver parágrafo se- 
guinte) pode possuir uma vu 12 saídas, o registador 
pode ter uma ou mais entradas. Mas como nos interessa 
uniformizar este circuito com o anterior, temos um con- 
versor paralelo-série identico ao descrito anteriormente, 
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dae modo a termos uma única entrada séria. Nestas con- 
dições, o registador de saída será idêntico ao anterior 
(agora o circuito integrado é constituído por 25 F.F. 
como se indica na Fig. n.º 7). 

A entrada do digito de controla (1) é feita no co- 
maço de cada sequência em que o CP faz o Set auto- 
mático para zero do 1.º F.F. e, o aparecimento do pri- 
meiro «dado» ao 2.º impulso leva o Set à posição nor- 
mal. 

No fim de cada sequência, este bit de controle 
não só nos vai abrir as saídas dos F. F, como também 
fazer o Set para zero do 1.º F. F. preparando-se o regis- 
tador, que está «limpo» para receber um novo fluxo 
de informação. 

Como a corrente de saida dos núcleos de memó- 
ria é bastante grande (3000-500 mÃ) comparada com a 
corrente admissível por estes circuitos integrados (da 
ordem do mA), à entrada do conversor paralelo-séria 
encontram-se resistências limitadoras de corrente de 
modo a baixar o valor de corrente de entrada. Por 
axemplo, se a corrente de entrada passa da ordem dos 
3-4 mA em cada entrada, teriamos um sistema cons- 
tituido por resistências na proporção de 1-100 associa- 
dos em paralalo: 


400 mA 100 KA 
€ ssa BUFFER 
- R7 — 
| Ps 
| ” À |2 
1 
R1 IK AL 
R1 
o EE 
R1 + R2 


1 
H=4x107T!— 4 mA ied. 
101 


O sistema de comando para a abertura de saídas 
desta unidade já vimos como era feito, 

Vamos agora indicar como se processa o sistema 
da comando. 

No início de cada sequência, existe um START 
automático que faz o SET para zero do 1.º F. F. Assim, 
à chegada do 1.º bit de informação, o F.F. já está 
ocupado, pelo que, a informação aí armazenada é trans- 
ferida para o 2.º F.F. O próximo bit provoca um deslo- 
camento e assim sucessivamente, 

Quando a informação estiver complata, (12 ou 24 
bits, consoante se trate do registador de entrada ou 
saida) o bit «1» de controle transita para o 13.º F.F. 

As saídas de todos os F.F. vão dar a portas AND 
as quais só abrem quando aparece um impulso «1» 
no 13.º F. F. Nestas condições, como o bit de controle é 
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«li», o é armazenado ressa unidade, é esse mesmo bit 
que provoca a abertura das portas e a saída da infor- 
mação. 

Mas, como não nos interessa manter as «portas» 
sempre abertas, mas sim fechadas, quando aparece 
um novo fluxo de informação, um «1» é armazenado 
no 1.º F.F. para começo de uma nova sequência € 
todos os outros F. F. se encontram «limpos», 

Nestas condições, a saída deste 13.º F.F. vai 
actuar sobre um circuito monostável, cuja saida é «il» 
(devido a um inversor) o qual vai fazer o RESET para 
zero de todos os F.F., excepto do 1.º. Neste, o SET 
é feito para zero num instante At (ou mais) em relação 
ao fluxo de informação de saída de modo a que quando 
os F.F. forem a zero já este fluxo tenha passado para 
a unidade seguinte, 

Esta condição também poderia ter sido assegu- 
rada por um circuito DELAY regulável. 

Se não quisesse utilizar este tipo de controla em 
vez de utilizar 13 F.F., poder-se-ia só utilizar 12 F.F,., 
em que o comando das portas seria efectuado por um 
contador-por-doze (ao fim de doze impulsos dava uma 
saida alta) que abriria quando a informação estivesse 
já armazenada nos 12 F.F. 

Se por acaso não quiséssemos utilizar nenhum 
destes sistemas, poderiamos então utilizar um convear- 
sor paralelo-séria (fig. 8). 


semi a . : 
: | eee 
; a Lo dgtb 
+ & 
[] n | 
| Lito 28 | | 
& 
| mis SR a | 
izisss se. Placa a a a ==. 4 
A | 
Ego mtrio id 


| SMpA Do 437 EF 


Saída do monostável que liga o RESET do 2.º até 
ao 13.º F.F. e ao SET do 1.º F.F. 


Este circuito é constituido por um registador de 
deslocamento série-paralelo em que a entrada-séris é 
designada por |, e a saída-paralelo possui 12 F.F. de- 
signados por: 
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LL. k, JH, GF, E D, E B, A, por um gerador 
de impulsos 'O.,, com a mesma frequência do conversor 
paralelo-séria já descrito, e por um conjunto de elemen- 
tos para o controle automático do conversor. O regis- 
tador de deslocamento é constituido por 13 elementos 
sendo os 12 primeiros destinados ao armazenamento 
da informação e o 13.º a funções de controle. O gerador 
de impulsos é constituido de um «SCHMITT TRIGGER». 
o conjunto de controle compreende um F. F. R. 5. básico 
designado por F. F., um trinco bi-estável designado por 
FF., dois circuitos AND, E, e E. e um inversor IN. 


Para a análise do conservador série-paralelo aci- 
ma descrito partiremos de um estado inicial caracta- 
rizados pelas condições 0, =O0O S=ÀAÁ= = 
= O e pela ausência da informação em |, (entrada- 
-séria) o que implica S= O em FF,; a condição So 
=0=>R=õOnomesmoFF,. À condição Q,= O blo- 
queia a transferência para CP através do circuito lógico 
«AND» E dos impulsos gerados por um O,. 


Como L= O, a mesma condição faz com qua o 
valor lógico O seja constantemente transferido para a 
saida O. de FF. pela a acção contante em T dos impulsos 
gerados em O, . Então O: = O e como R, = O: é R = 
O, mantendo-se, neste estado inicial, o regitador com- 
pletamente «limpos». 


Consideramos o seguinte fluxo de informação: 


010 100 000 000 


a qual é precedida do bit 1 de controle, Logo que este 
bit chegua à entrada do conversor, como R = 1, a en- 
trada de S de FF, assume o valor lógico 1. Então a con- 
dição S = faz mudar de estado este F. F. o que se traduz 
na mudança de O,= O para O, = 1. 

Com L= OQ, = 1, passa agora a ser constantemente 
transferido para O. = 1 condição que prepara o regista- 
dor para entrada de informação. Ao mesmo tempo, a 
mesma condição OQ, = 1 desbloqueia o circuito lógico F.. 
Nestas condições, começam ocorrer impulsos em CP 
que fazem entrar, passo a passo a informação série 
precedida do bit 1 de controle. Ao terminar 13.º im- 
pulso, a informação apresenta-se na forma paralela nas 
saidas La À, eo bit 1 de controle em So - Admite-se 
ainda que estas saídas estão ligadas às entradas de um 
registador de transferência de memória RTM constituido 
por F. F. emaster/slaver e que a saída S, está ligada à 
entrada CP do mesmo RTM. Estão, a condição S,= 1 
faz transferir para as secções «master» dos F.F, do 
RTM. os valores de L— a — À, ficando retida a infor- 
mação, neste RTM, que será, no 1.º caso, (registador 
de entrada) a matriz distribuidora e os núcleos de me- 
mória e, no segundo caso o conversador D/A, 


Ao mesmo tempo qua faz esta transferência, como 
R = So a condição S;= 1 faz com que R=1 em FF, 
condição que acarreta S = 0 e O, = O. Esta condi- 
ção O, = O, bloqueia o fornecimento do 15.º impulso 
a dos subsequentes a CP, e, como L= Q,, a ocorrência 
do 15.º impulso em | faz transferir o valor lógico «O» 
para Os. 
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Com O.= 0, e Rs = O, vem R, = O, que «limpa» 
o registador. Esta «limpeza» inclui a mudança de So = 1 
para So = O o que faz transferir completamente a infor- 
mação para as próximas unidades. 

Entretanto como ocorre na análise precedente, o 
conversor fica de novo preparado para receber um novo 
fluxo de informação. Qualquer dos métodos descritos 
anteriormente realiza plenamente o circuito lógico dese- 
jado, pelo que a escolha deverá recair no mais preciso. 
Contudo, esta característica não é reconhecida facil- 
mente pois, os sistemas apresentados são bastante 
semelhantes. 


2.3.5 — INTERRUPTORES MULTIPLAS 


SAÍDAS 


COM 


Uma parte muito importante do nosso computador, 
a matriz de comutação, vai ser estudada neste mo- 
mento. 


X Y Ci S Co 
o ol O O 
DO 1 t O 
Oo 110 1 0 
oO 1 1 O 71 
1 0 0 1 0 
TO o 1 
1 1 0 O 1 
1 MM A 1 4 


Desta tábua tiramos: 
S=XYCI+XYCI+XYCi+XY Ci 
Co=XYCi+XYCI+XYCI+XYcCi 

em que X, Y, Ci são respectivamente os bits adendos é 


o bit de transporte. Podemos escrever S e Co em fun- 
ção dos termos P, 


S=ALR+R+P, 
| Co=P.+P+P+P, 


S = (OxPo) + (1xP,) + (1xP,) + (0xP.) + 
+ (1xP,) + (0xP;) + (OxP,) + 


+ (1xP,) 
=| 
Co = (0xP,) + (0xP,) + (OxP.) + (0xP,) +- 


+ (OxP,) + (1xP;) + (1xP;) + 
| + (1xP,) 


Estas duas últimas equações podem escrever-se 
na forma de uma matriz Booleana: 


s e] Ps 


00010111. P, 
Co Ps 
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ou na notação: 
DO] = [B] . [P] 


em que O, B, P são matrizes. Contudo podemos dar 
outro aspecto às funções S e Co 


o =— Do 3, Ss. Sa 
—— pe 
Co = S 5. Se. 5 


10 +S) (0+S) (048) (1458) 
East 
Co=(0+S) (0+8S|) (0+S.) 


1+S) (M0+S) (1+S)(1+S) (1+ 5) 


ou na forma de matriz: 


87] 01191001 Sa 


90010111 5: 
“Co Sa 


ou na notação: 
[0] — [B] [5] 


Como se vê em qualquer destes dois exemplos temos 
uma matriz de entrada, uma matriz de transferência 
a uma matriz de saida, sendo as saidas em número mul- 
tiplo das entradas. Servimo-nos desta propriedade para 
conseguir várias saidas à custa de poucas entradas, 
utilizando uma matriz de interruptores. Dos vários tipos 
de matrizes de interruptores conhecidos vamos usar as 
matrizes de diodos. O número de entradas necessárias 
para o nosso circuito são doze, o que dá 4096 saidas 
(para n entradas há 2n saidas). Sejam as entradas 


a, & ..., dj: & 05 respectivos conjugados. Para este 
conjunto de entradas teremos as saidas fi, fi, ..., Tue, 
que podem ser escritas por função Booleanas. 
f — a a . a. - mi. a 
1 j* Fã 3 “12 
E. = a “ a. a RR = 
Ê 1 É 5 12 
Ê = a O a “ a. .... «ã . 
3 1 ê 3 11 Je 
E » as = à Bol. a, 1 ves à A 
4096 1 Ê 11 le 
F = a . d . d.. “o. ” 1 A a Ba 
EildZ0 1 i- + 1 Í 1 E 


150 


ou sob a forma de matriz: 


r | | ca 
| : 1 | = | | L | Ê O à à =. Ú O [Fm 
d + P | | 
Ls ho Hj tl [] E E = O O P? 
[3 UU U | DO 0 O FP 
e r] Ú ( R| 1 WU à () U Pd 
E 
ams Dus cd QD. a .. o O FP: 
n * * a * 
E * * “ * * * 
m = = LL] 
ei 8) M) U t Le » a a ” 1 ( Pá Õ “q + 
! al) +93 
GE RU SE | 
| O O BO Disso O 11 /P4096) 
“4039 E — A 


Nesta matriz Booleana aparece «1» em cada linha 
e todos os «1» pertencem à diagonal da matriz. Como 
se vê só existe uma saida capaz por cada linha, donde 
se conclui que a matriz de diodos funciona como uma 
porta AND, O número de diodos necessários para rea- 
lizar esta matriz é de 49.152 (para n entradas nx2 1 
diodos). Contudo podemos simplificar este número uti- 
lizando matrizes de árvores duplas, Nestas condições 
as equações de saida tomam o seguinte aspecto: 


Um fio a atravessar o núcleo de ferrite com uma 
corrente Im teriamos 2 fios a atravessarem esse núcleo 
cada um deles percorridos por uma corrente Im/2 (ver 
figura n.º 9b). 


Quer esta solução, quer a anterior, vão dar saídas 
em excesso, pelo que teremos de ter um relógio sincro- 
nizado com os elementos fornecidos à memória, de 
modo, a que haja interrupção logo que se atinja o nú- 
mero de bits que é possivel armazenar na memória 
(paragem do relógio). 

Para a primeira solução, utilizam-se os 4.000 nú- 
cleos da memória, o que corresponde a 333 séries 
da entradas de 12 bits que ao fim destas 333 suces- 
soes não há mais fornecimento de dados ao registador 
ou desliga-se este automaticamente. 

Para a segunda solução vai-se utilizar uma matriz 
64x64, mas como a matriz da memória a utilizar é 
50x48 e outra 16x16, ficando 144) núcleos livres. 
Esta solução só permite 213 séries de entradas de 12 
bits. 

A 1.º solução admite 32 séries de 25 minutos e a 
segundo 16 séries de 25 minutos, pois por minuto trans- 
mite-se uma série de 12 bits. 


2.3.6. — Matriz selectora 


Esta matriz possui dois grupos de terminais de 
antrada e um de saida. 


TÉCNICA 423 


MATRIZ 
SELECTORA 


= Tl 
º entradas 


1 saída 


n entradas 


Um dos grupos consiste em 2" terminais de en- 
trada que estão directamente ligados às «cabeças» lai- 
toras de memória. Os outros n terminais de entrada ser- 
vem para selecionar o canal desejado. O terminal de 
saida de 24 bits. Esta matriz selectora tem como fim 
transportar para o registador de saídas dos dados arma- 
zenados na memória pela mesma ordem em que estes 
aí se encontravam, isto é, lã os dados que se encontram 
armazenados na memória que previamente foram selec- 
cionados pelos canais selectores. Comparando o tipo de 
matriz que temos à saída da memória (para ler os da- 
dos) com a que tinhamos à entrada (escrever os dados) 
notamos que a segunda é a inversa da primeira, com a 
única diferença dos dados agora sairem em série, por- 
menor que se quizer pode ser eliminado. Nestas condi- 
ções pode-se fazer a utilização dupla da matriz anterior, 
utilizando um interruptor que, primeiramente, coloque a 
matriz ligada na posição de escrita ql» e que depois 
comute para uma posição «0» de leitura, invertendo a 
polaridade da fonte da polarização da matriz. 


Assim se considerarmos a primeira solução, se a 
escrita se fizer pelos fios Y que atravessam a memória 
(ver figura 9a) a leitura far-se-ã pelos fios X ou vice- 
versa. Para a segunda solução far-se-á a leitura por Z 
a N. Contudo a melhor solução será a que utilize as mes- 
mas linhas quer para a entrada e saída, utilizando um 
interruptor único, que inverte a polaridade da fonte de 
polarização da matriz e a torna capaz de conduzir em 
sentido contrário, quando o núcleos da memória são 
excitados. O papel dos selectores de canais é precisa- 
mente determinar quais as colunas da matriz da meméó- 
ria escolhida para desexcitação dos respectivos nú- 
cleos. Asim,cada um dos 12 canais abrange 4 colunas 
(ver figura 10) as quais são seleccionadas por meio de 
uma porta de comando, que selecciona o canal a utilizar, 
operação esta que que se pode fazer manual ou automá- 
ticamente,por meio dum impulso de relógio com um At 
perfeitamenta sincronizado. 


Uma vez seleccionado o canal, existe um F.F. D 
ligado a duas das colunas com um At tal que permita 
a leitura dos 12 bits da primeira coluna no primeiro 
ciclo, e a leitura dos 12 bits da segunda, no ciclo 
seguinte, operação esta que se processa até à leitura 
de todos os bits das duas colunas mencionadas. Uma 
vez lidos todos os bits, existe no fim de cada sequência, 
um núcleo no qual existe um «1x armazenado e que ao 
ser excitado permite a abertura de uma outra porta que 
dá passagem ao impulso T para o trigger seguinte e 
provoca o STOP do primeiro trigger fazendo o SET para 
o zero deste. 
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À segunda séries de valores das terceiras e quartas 
colunas é lida por um processo idêntico, simplesmente 
agora o «il» armazenado no enúcleo de comando» val 
provocar o STOP do primeiro canal fazendo o RESET 
automático, ou o SET para zero e, o SET para «il» do 
2.º canal que por sua vez vai actuar sobre uma porta de 
comando que irá seleccionar as 5.º, 6.º, 7.º e 8.* colunas 
e assim sucessivamente, 


Por outro lado, ligado com o T de cada F. F., tere- 
mos um contador de 12 bits, que vai actuar sobre uma 
porta de comando à saída dos registadores que permite 
a abertura da mesma e a passagem simultânea dos 
dados da altura e azimute para os conversores analó- 
gicos respectivos (figura nr. 10). 

O núcleo de comando da 50,º coluna (1.º solução) 
ou 64.º coluna (2.º solução) vai actuar sobre o comando 
do relógio de modo a fazer STOP automático de todo 
o circuito e do comando da matriz de modo a colocá-la 
de novo na posição de «escrever» uma vez que a memó- 
ria esta pronta a receber novos dados (ver figura nr. 10). 


—— Bt 


p= 
h az b 
f j (é 


: E A d 
contados : q telas 


fo test A qua) 


Nota: — Pode também usar o selector de canais 
para a primeira solução do comutador a utilizar desde 
que a corrente à saída da matriz seja Im/2. Nestas 
condições a outra corrente Im/2 iria passar no selector 
de canais cuja utilização é dupla. 


2.3.7 — Unidade de memória 


A unidade de memória por nós utilizada é consti- 
tuida por quatro planos úteis que possuem 4.000 núcleos 
dispostos numa matriz 50x80. Esta memória possue 
duas linhas de percepção que atravessam os 2.000 
núcleos de uma matriz de 50 x 40. Duas linhas inibi- 
doras paralelas às linhas de percepção atravessam 
também os núcleos da memória e estão ligados em 
sério em forma de W. Além destas linhas temos de 
considerar as linhas leitoras e de escrita que vão ler e 
armazenar respectivamente os dados. 


Vejamos como é constituido cada um dos núcleos 
da matriz. 
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2.3.7.1 — Caracteristicas do nucleo magnético 


O núclso magnético de ferrita possui uma curva 
de histeres e como se indica na figura 11 - a, que nos 
dá a grandeza do campo magnético, H, e do campo de 
indução , B, do material, Quando a intensidade do cam- 
po magnético é núla, exista uma densidade de fluxo 
remanescenta Br chamado fluxo de indução residual. 
Portanto, na memória, existe dois estados binários: O 
fluxo remanescente positivo, + Br, que pode ser desi- 
gnado pelo estado binário «1» ou SET e o fluxo rema- 
nescente negativo — Br, que pode ser designado pelo 
estado binário «O» ou RESET. 

Para variar o fluxo remanescente deve-se aplicar 
um campo magnético de intensidade Hc chamado «cam- 
po de força compulsiva». Um núcleo pode ser comutado 
de um estado para o outro por aplicação de um impulso 
de corrente atravéz de um dos fios que o atravessa. 
Durante o 1.º meio ciclo da curva de histerese (curva 
mais a cheio da figura n.º 11 -b) o núcleo passa de «0» 
para «ly por aplicação de uma corrente Im. No meio 
ciclo seguinte comuta de «ly para «0». Nestas condi- 
ções, a transição do estado da memória, é comandado 
pelo impulso de corrente Im que nos conduz a um valor 
da Hc máximo, de modo que a curva de histerese se 
aproxime o mais possível da forma rectângular, o que 
diminui a constante de tempo de comutação que nos é 
dada por: 


G = para H > 2 Hc 


Ho —- valor do campo magnético no instante—= O, 


o qual é determinado por extrapolação quando 1 —> 


H— é valor do campo magnético no Instante 
considerado 


Sw — coeficiente de comutação que nos é dado 
pela inclinação da recta 


Sw = Se + Sr em que 
58 — representa a perca de corrente que reflui 


or — representa a relaxação do circuito 
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Numa primeira aproximação podemos dizer que a 
perca de histerese + a de corrente que reflui no circuto 
+ a relaxação do circuito dá: 


4 Im 
4 Ho Im + Sw — 


Bm 
em que Im é dado por Ea 


Como se pode deduzir da exposição interessa-nos 
saber qual a corrente de magnetização e tensão a que 
se deve submeter um núcleo quer durante o armazena- 
mento, quer durante a leitura, Nestas condições, inte- 
ressa-nos saber qual a relação entre a tensão de entrada 
(armazenamento) e de saida (leitura) do núcleo, as quais 
vão determinar as polaridades dos diodos da matriz 
comutação fazendo com que estas conduzam, ou não. 


Supondo que a tensão de alimentação do circuito 
é de V e que a resistência dos diodos é R1 à entrada 
a R2 à saida temos que: 


Para os núcleos de ferrite em causa a corrente ce 
magnetização é de cerca de 500 a 600 mili-Ampéres 
a é 1 us coma = kI%mvV,káÁssim tem-se que e, 
é a tensão induzida aos núcleos durante o tempo de 
comutação. À resistência dos diodos é da ordem de 1 
a 0,5 ohms. 


Das condições anteriores podemos tirar que: 


em que(Z! e(7) m são os fluxos residual e máximo res- 
pectivamente e vêm em maxwells. Seguidamente ilustra- 
mos um caso de armazenamento de dados na memória 
em função das expressões anteriores. 
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Tendo em conta os valores de e; e er temos que 
ArtQ)m é 3 maxwells, e atensão V de excitação da 
matriz (tendo em conta a baixa resistência da matriz e 
do próprio núcleo) é de 12 a 15 volts (R, — 20 a 
23 ohms). 

Como para as operações de leitura e escrita se 
utiliza a mesma matriz conclui-se que e:= — e, isto 
é, basta polarizar a matriz em sentido inverso e fazer 
com que i. = i;/2. Normalmente quer na escrita quer 
na leitura poder-se-iam duas correntes iguais a metade 
do valor máximo que passariam através dos núcleos 
magnéticos, impondo excitação dos mesmos. 

Para o nosso caso a operação de leitura procussa-se 
segundo um esquema semelhante ao indicado, 


2.3.7.2 — Armazenando da informação nos núcleos 
magnéticos 


O armazenamento dos núcleos de ferrito é feito 
segundo a direcção de magnetização destes. Assim, se 
o fluxo magnético tiver o sentido 1 (figura 11) dizemos 


que o núcleo armazenou um 1; se tiver o sentido con- 
trário dizemos que armazenou um O, 


Pela mesma figura vê-se que o ciclo de histerese 
da ferrite é grande, e, se fluir no núcleo um impulso 
de corrente de amplitude + If (designada por ful 1) 
então é induzido um campo magnético da amplitude 
+ Hm. A densidade de fluxo no núcleo sobe a + Bm e 
permanece no valor remanescente + Br quando o im- 
pulso de corrente a o campo magnético são removidos. 
Se, no entanto, fluir no núcleo um impulso de corrente 
(— If) então o campo induzido será — Hm e a densidade 
de fluxo mudará de + Br para — Bm e permanece em 
— Br quando o campo é removido (figura 11). Se se 
aplicar de novo um campo — Hm, a densidade do fluxo 
passa simplesmente de — a — Bm. 

A corrente + If, é designada por «correnta de 
escritar e levadaa núcleo ao estado 1 (+ Br); a cor- 
rente — If é designada por «corrente de leitura» e leva 
o núcleo ao estado O (— Br). Se uma «corrente de 
leitura», — If é aplicada a um núcleo que está no 
estado 1, a designação do fluxo muda de + Br para 
— Br; se o núcleo estiver no estado O, a mudança faz- 
-se de — Br para — Em e em seguida para — Br. 

Assim, o estado de um núcleo, isto é, a polaridade 
da magnetização permanente pode ser determinada apli- 
cando um impulso «leitory a medindo a amplitude da 
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f. e. m. induzida quando a densidade do fluxo muda, 
Observando o ciclo de histerese, uma meia corrente 
+. |f/2 não comuta um núcleo dum estado para o outro 
pois a densidade de fluxo varia de 1 B, e volta a 4 Br. 
A aplicação simultânea de duas meias correntes no 
mesmo sentido, terá o mesmo efeito que uma corrente If 
nesse sentido, o que pode ocasionar a comutação de 
um núcleo, 

Usam-se estas propriedades de «meias correntes 
coincidentes» para a selecção de um núcleo específico 
da memória. 

Vemos assim que a corrente leitura — If é destru- 
tiva pois o núcleo é conduzido ao estado O. Se a infor- 


mação armazenada num dado núcleo deve ser retirada 
depois da leitura dum 1, é necessário retornar a rescre- 
vê-lo depois desta leitura, 

A memória construída por este processo é mais 
económica do que utilizando registadores. 

Assim para seleccionar um núcleo, o princípio das 
meias correntes pode ser adoptado, baseando-se no 
facto de que um núcigo não pode ser comutado por uma 
meia corrente. Dois fios, um de escrita-X-horizontal e 
um de leitura — Y — vertical são entrelaçados como se 
mostra na Fig. 12. Além destes fios, há o fio S (Sense 
Wire) através do qual é medido o sentido da amplitude 
da tensão induzida, que se encontra enfiado nos nú- 
claeos em diagonal, Há ainda outro fio, Z, fio inibidor 
que se encontra nos núcleos como indica a Fig. 12. 
Este fio é usado durante a operação de escrita, em par- 
ticular quando vários núcleos da matriz são usados em 
paralelo. À informação pode ser armazenada na memória 
por meio de duas meias correntes (1/2) que fluem atra- 
vês dos fios X e Y independentemente do estado do 
núclao. Contudo, se se pretenda armazenar um O basta 


O 


fazer fluir pelos fios uma corrante negativa. 

Pela Fig. 12 vê-se que cada elemento da matriz 
é constituido por quatro núcleos de ferrite que arma- 
zenam um único bit octal, 


2.3.7.3 — Operação de leitura 


Esta operação é feita através da matriz selectora 
segundo o esquema temporal indicado anteriormente. 
Primeiro, tal como se indica na Fig. 13-a, uma corrente 
leitora é conduzida através dos fios Xm, Yn para a ex- 
citação dos núcleos correspondentes de modo a obter-se 
a informação que estes continham. Quando se lê um O 
as tensões induzidas são aplicadas a portas AND co- 
mandadas por um impulso. Nestas condições as portas 
ficam bloqueadas indo para o registador a informação 0, 

Normalmente a fonte de tensão e a tensão indu- 
zidas estão desfasadas no tempo, fazendo com que haja 
um atraso na comutação do núclso o que origina uma 
desfazagem nas correntes leitoras. 


2.3.7.4 — Operação de escrita 


Como se disse anteriormente o armazenamento dos 
bits na memória faz-se por meio da matriz comutadora 
segundo qualquer dos dois processos indicados ante- 
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REA GOVACE oa O O O 


4 +30 lv) aa, 


4 


READ sina CT) q 


4 Bold) 


Y READ T, e | | | | 


4509 


reta hã o 


riormente. Para qualquer um destes esquemas o arma- 
zeramento é feito directament por uma corrente Im. 
Em vez de se utilizar um dos processos anteriores 
poderiamos utilizar um terceiro, que consistiria em fa- 
zer-se uso simultâneo da matriz comutadora e do cir- 
cuito (agulheira) dos canais. Nestas condições o reló- 
gio do registador de entrada deve estar sincronizado com 
o dos F.F. dos canais. 

Em seguida indica-se o esquema temporal de es- 
crita que é válido para qualquer dos três processos 
descritos 


processos descritos. 


Gini 1 
dt FO a 
= els ml 
A GAME atgina TT | Es preta tos 
Ko San so ol 


a O [mem A 4 


Ke DE E o O 
e sw asse ii DO 
E sidk T | | | 

IL T 


Na Fig. 13-b indica-se um esquema de armazena. 
mento de 4 bits. Na Fig. 14 indica-se a esquematiza- 
ção para os dois modos de operar os núcleos de me- 
mória, 


RATE uMsTÊ 


2.4 — Circuito subtrator 
Este circuito vai permitir com que as sequências 


de saída sejam tais que a acção sobre a embraiagem 
só se faça sentir pela diferença entre a sequência pre- 


tod 


sente e a antecedente, fazendo com que a antena se 
desloque num ou noutro sentido. O circuito que vai rea- 
lizar tal é o seguinte; 


saida | 
para O buftel 


R, R: R: — Registadores R— 13 FF 
C' — Circuitos complementadores R.e— 13 FF 
F. A. — FULL-ADDERS R.— 14 FF 


D. E, — Circuitos DELAY 
REGISTADOR 


Data lr | DE E 
$| EF, FF: Ffy | 


COMPLEMENTADOR(C,") 


pos pe 
FF F Fa 


No caso de termos a sequência 0110001010 pro- 
cedida ou seguida na sequência 0010110111 equiva- 
lente aos nrs. (612), e (267), à saída do 1.º comple- 
mentador teremos: 
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a) (612) + (— 267), 


e ms e 


A, B 


O -> positivo 


1 — negativo 
Portanto teremos o valor complementado B 

(0010110111) = 1 101 091 007 

Em C, teremos: 


o 112 001 010 
+ 1 101 001 091 
xX'0 011 010 011 
| 


cespreza-sa 


- (323).c.e.d. 


Como o «indice-sinalh é zero, o resultado é posl- 
tivo pelo que a porta 1 permanece fechada e a porta 2 
aberta. Ão 13.º C.P. deslocamos a informação para 5, 
[sistema análogo aos indicados para os registadores) 
sendo este o valor que se encontra no registador de 
saida. 


b) (267), |— 612), 


.— 


A 8 
(0 110 091 010) = 1001 110 191 


E em C,teremos: 


O OO 110 111 
1/09 110 100 


1 100 101 100 


O sinal é negativo, logo fica aberto e 2 fechado, 
pelo que teremos 5, = O constantemente, e 5. igual a: 


O 011 010 011 i.e. (323), com sinal menos 


Como se vê, C: será um simples inversor, portanto, 
muito menos complicado que C, cujo circuito se repre- 
senta em (C). 

Para o nosso caso, existe à saida da memória um 
biestável com um tempo de comutação, sincronizado 
com o monostável de comando de total modo que os 
primeiros 12 bits (altura) serão «encaminhados» para 
um subtractor do tipo descrito e a segunda sequência 
para outro subtractor em que não existe a unidade C'. 
pois a diferença entre o dado presente e o antecedente 
nunca é negativo. Nestas condições prescindimos tam- 
bém das portas 1,2 e 3. 


O circuito funciona do seguinte modo: 


Ao seleccinarmos o «canal de memória» accio- 
namos também o comando do biestável o qual ao fim 
do tempo 7, equivalente à 1.º sequência de 12 bits 
comuta para a posição dois durante um tempo também 7 
e assim sucessivamente. Ássim, a 1.º sequênca de 12 
bits passa ao registador R, e R: encontra-se apagado. 
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Nestas condições, o valor de R, soma-se ao de R: 
e dá-nos zero, ao mesmo tempo que a informação de 
R, é armazenada em R. depois de sofrer uma «comple- 
mentação a dois». 

Assim, à passagem do 1.º bit CP está em zero, i.e. 
no início da contagem, e o circuito complementador 
“prepara-ser para começar a receber a informação e 
complementá-la, para depois registála em R., a qual 
vai ser subtraida através do F, A. à próxima sequência, 

No caso do azimute a informação recolhida à saída 
do F. A, é armazenada em R; que neste caso não é 
mais que o nosso BUFFER. No caso da altura, é neces- 
sário um 2.º circuito «comparador» como se indica, e é 
saída desse comparador que vai ser armazenada no 
BUFFER correspondente. 

As saidas destes BUFFERS são automatizados de 
modo a que a informação passe simultansamente do 
BUFFER para o C. D./A./T. 


2.5 — Conversor digital-analógico temporal 


Como já se disse anteriormente, uma vez lidos os 
bits correspondentes à altura e azimute dos respectivos 
núcleos de memória, surge um comando que actua so- 
bre as portas à saída do registador abrindo-as, de modo 
a que a informação transite para o conversor D/A.T. 

Nestas condições, vamos tentar construir um cir- 
cuito que converta o sinal do registador num dado nível 
a fim de se obter um tempo de actuação concordante so- 
bra a «embraiagem» de comando do motor, que acciona 
a antena helicoidal, 


O nível de tensão correspondente ao valor lógico 
«ly à saida dos registadores é da + 5 €— + 10 MV. 
Baseando-nos neste valor, e no nível de actuação da 
embraiagem de comando (— 5V) e no tempo necessá- 
rio de actuação (da ordem dos segundos), vamos pro- 
jectar num sistema que satisfaça a todos estes reque- 
sitos. 

Primeiramente, vamos estudar um sistema «descri- 
minador» das tensões, consoante a ordem de entrada 
dos bits, do mais significativo ao menos significativo. 


2.5.1 — Conversos dgital analógico 


O conversor vai ser formado à custa de resistân- 
cias de precisão. Assim, vamos explicar o funcionamento 
para um par de quatro entrada considerando primeira- 
mente o efeito da existência dum nível alto (os outros 
são baixos). Depois, recorrendo ao princípio da sobre- 
posição, apontar-se-á uma solução final para todo o 
sistema de entrada. 

Nesta Fig., pusémos propositadamente interruptores 
nos nós 1, 2, 3, e 4 para significar a presença ou 
ausência de impulso. 

Vamos começar por determinar a resistência de 
saida (Ro) do circuito de posição de OFF, começando 
por analizar o n.º 7. 
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——s mae (6) 


214% 
- [Ea 


| - 


Neste N.º tem-se: 


1 


>>> = R 
1/8; + 1/Rs 
1 1 
NÓ 8 ESP E = R 
1/R: + 1/R, + R 2/2R 
1 
NÓ 5 PQ e E 


1/R, + 1/R:+ BR) 


O valor destas resistências R vai ser função da 
corrente de condução, 

Consideremos agora a existência de um «ly só na 
porta 1, |. e. a existência da ordem mais significativa. 
Nestas condições, poderemos transformar o nosso 
esquema no seguinte (equivalenta): 


VoA 


RL 


Nestas “ondições, se considerarmos R, >>R 
tem-se: 


2R 1 
Vo, = a V —Y 
da Nona A 2 "R 
a, valor de R; será 
1 RL 
Vo mm PS 
Aê Ria +R 
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4) 
: va 


[ESQUEMA EQUIVALENTE 
ÃO ANTERIOR) 


Consideremos agora só a existência de tensão no 
nó 2 (próximo bit mais significativo): 
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Relacionando-sa esta tensão com Voa tem-se: 


para um R, qualquer tem-se: 


Supondo agora a existência de tensão só em (3) 
tem-se:: 


v2 
E.» 
21 o. 
“4 º|2+22., io Va 
21 
6 
pm: É 
5 6 
Va = MV, Dec, R 
6 R+R 32 
2 R 2 6 1 
Vo, = 'l —) MV | === a o 
A 2 R+R ' 3) g A 
1 Rt 
Vo SEO Va (——) 
à E a 


Assim, se quisermos determinar qual a solução para 
as três primeiras entradas iguais a «1» bastava-nos 
aplicar o princípio da sobreposição aos valores deter- 
minados: 


1 1 1 Ri 
E 
7 Re 
“8 'StprR 
7 
se R|, >2>2R5> Voy Va 


Assim, se aumentassemos o número de entradas 
para quatro, opróximo bits corresponderiaa 1 Vp de 


16 
tansão à saída e assim sucessivamente. 
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Nestas condições, para qualquer número de entra- 
das teremos um valor de tensão máxima dada por: 


Va 
Vo= (21-11) — 
2" 


Assim, podemos escrever uma equação geral que 
satisfaça todos os valores de entrada i.e.: 


1 1 
Vo.=(— Db, + -— Ds; + + TO 
da 2 4 8 
1 RL 
+—= DI) X V ERES 


Em que D representava o valor lógico da entrada, 
se é qi» ou zero. 

As vantagens em utilizarmos resistências de Valor 
R e 2R saltam agora à vista, pois conduzem-nos a uma 
expressão relativamente fácil, em que a resistência de 
saída do circuito (Ro) é constantemente igual a R. 


Para o nosso caso, n = 12. 


1 1 1 
nm Do si o si ce [d 
Sr as as 1 ya a) 
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Carta diois di Casuparadar, e Stuatah a 


fLE SET Do 


(iam ' Crua TO 


R aguivadadt: Rpg * Rat Ra + RE 


Eur EuTfasTARE 


entrada positiva corta a — 4 V. 


Quando o integrador atingir esta tensão, O compa- 
rador de saída comuta para o estado tl» e o transistor 
O, não só reproduz o impulso de entrada com a dura- 
ção do tempo de comutação, como também faz o 
RESET da rampa do circuito integrador (gerador da 
rampa); define-se assim o começo da amostragem da 
rampa negativa e o seu fim, assegurando uma tensão 
de O V à entrada (+) do comprador até ao fim do 
intervalo de RESET. 

O circuito, fica livre para começar um novo ciclo, 
Porque o intervalo do RESET é muito rápido, quando 
comparado com o tempo da rampa, a relação de repe- 
tição dos impulsos de saida é definida quase que intei- 
ramente pelo período da rampa do integrador (a qual 
só exista quando houver o sinal de entrada) deixando 
as partes de manutenção para as partes restantes do 
circuito. 

No começo da contagem dos intervalos de tempo 
deve-se ter: 


1— À saida do integrador OV 

2 — O comparador na posição OFF; 

3 — À entrada (+) do comparador a — 4 V,, de- 
vido às resistências R, e R. (divisoras de ten- 
são). 


Nestas condições, a rampa gerada à saida tem uma 
inclinação negativa e, é traduzida pela seguinte equa- 
ção: 


à Vout Vin 
O (1) 
a T R, C, 

Como a tensão de saida é a correspondente ao dis- 
paro do comparador, i.s. quando ele passa da posição 
OFF a ON i.e. a entrada (—) do comparador passa a 
(+), e isso dá-se a — 4V, podemos pôr: 


4 R1 C1 


Tn a 2. 
Vin ua 
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TEMPO DE RESET 
No começo tem-se: 


1 — À saída do integrador a — 4V 
2-— À saida do comparador de O, são ON. 


A primeira limitação na relação de descarga é 
feita por R.. Por outro lado, R, controla a corrente do 
colector no momento de comutar a OFF (no fim do in- 
tervalo de RESET) de modo a que o comparador comu- 
tará a OFF precisamente em OV, 


Assim, quando o comparador de saída comuta para 
(1) a corrente do emissor de Q,, determina no compa- 
rador uma entrada de tensão de referência (+) para 
comutar de — 4WV a aproximadamente OV. Ao mesmo 
tempo a corrente colectora de Q, vai-se somar ao borne 
(|—) do integrador, originando uma rampa estreita de 
saída neste gerador em RESET, para um dado tempo 
descrito por uma equação análoga a (2) em que R, 
é substituido por R*. Rf representa a soma de todas 
as resistências limitadoras de corrente na descarga do 
condensador integrador. 


R$=R 


1.º andar, saída T Ro PR 


= 250 + 170 + 25 +1300 
= 1685 9 


Fazendo Vin= 15 V e R6=R* em (2) tem-se: 


: 4RºC 1685 X 10-5 
Ta = —— =4z————— 
V — (fonte) 15 


— 449,667 5 


Atendendo a que, grande parte da limitação se 
deve a R, o 1.º, 2.º, 3.º termo de 3, são divididos por 
R. durante o RESET. 


Assim, para um Vin= 5V=> f= 12,5 Hz5T,= 
=4x10-2s 


Tº 


erre 
Tait Ta 


O erro do RESET (para este caso) = 


A este erro, podemos juntar um erro + 1% da 
fonte de alimentação e 0,5 % do tempo de comutação. 
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pr 


Volta 


ias; 


Para o nosso problema, as tensões de excitação vão 
variar entre 10 V. a mV. (Vin), pois a saida do conver- 
sor digital analógico vai ligar a entrada do conversor 
analógico temporal, 

Como se vê, quanto menor for o nivel de entrada, 
maior será o tempo de comutação, e, como queremos 
que o maior tempo corresponda ao bit mais significa- 
tivo, a ordem de entrada dos bits no conversor D/A 
deverá ser inversa da proposta. 

Por outro lado vemos que os termos de variação 
para excitação da embriagem do motor, variam entre 
um máximo de: 40 segundos e um minimo de 0,04 
segundos. Se por acaso se pretender aumentar o tempo 
de excitação deveremos colocar um monostável com um 
tempo de excitação considerando o mínimo (960 mili- 
segundos por exemplo) accionado pela saída do con- 
versor (um impulso, o qual se adiciona a este impulso 
por meio de uma porta «ou»). 

Para este caso, o fim do impulso deve sempre 
coincidir com o início do monostável. Como se viu, a 
impedância de entrada do conversor A/T é bastante 
elevada pelo que R (do D/A) é >> R e portanto 
é válida a aproximação feita durante o estudo deste 
conversor. 


E RR 
MONGESTANEL 


SALDA DO 
EXneiTO Dicyi dell 
ANALOGA vo 

TEM Toa L. 


mm 
| —— = 


Saida para 
o motor 


| MONO-ESTÁVEL 
1 (960 ns) 
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Também poderiamos arranjar o conversor D/A 
com uma precisão superior à indicada desde que 
se utilizasse também um amplificador direcional da 
ganho bastante elevado. Nestas condições teremos o 
asquema que seguidamente apresentamos: 


VY, = nivel que indica o sinal e pode ser O. 


VM. e V. = V, = Vo, — Va iai =Va 


Da equação tiramos: 


Vo V1 V. Va 
ma (SS 4 O Fu É ETs 
R 2R 4R 2,R 
1 1 1 1 
Gs (== 21== 4 + )V 
2 + 2, 2, x 


Nota — R é ordem das dezenas de Kº). 
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OBJECTIVOS E PROGRAMA 


Exposição Internacional de metal e seus derivados 


Foire Internacional de Machines et Materiaux d'Embales 


Congresso Internacional sobre Srtuctural Questions of Under- 
ground Railwais 

Secretariado: 

Kospek distudamanvyi 

Egvesúlet, 1655, 

Budapest 

Kossutk Lajos ter 6-B 

Hungary 


Exposição de Materiais Térmicos da Climatização e da Anti- 
poluição 


Colloque sur l'emplor de raymnents de haute intensite por 
la traitements des déchets 


Feira Internacional da Energia e da Petroquimica 


Feira Internacional de Lisboa 


Exposição Cientifica e Técnica 


Feira Internacional de Poznan 


Salon de L'éctronique Industrielle e de L'industrie életroté- 
chnique 


| 


| 
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Feiras, Conferências, Exposições e Congressos 


LOCAL, DATA, ORGANIZAÇÕES 
E INFORMAÇÕES 


Londres 
Abril - 75 


Dússeldort 9 a 15 — Maio - 75 
Interpack 75 


Hungray 


Scientific soc for transports 12 a 25 Maio 


Bulatou furea 


Lyon 


29 de Maio a 2 de Junho - 


Munich 17-21 Mais - 75 
Rep, Féd. d'Allemagne 


Teerão 
6a 11 Maio-75 


Lisboa 
11 a 22 de Junho - 75 


Grécia - Salonica 
6 a 12 de Junho - 75 


Polónia - Poznan 
8 a 17 de Junho - 75 


Bále - Suiça 
9 a 13 de Setembro - 75 
INELTEC 
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Equipamentos 


por ondas dirigidas 


destinados a transmissão dum canal de música ou de 
várias vias telefónicas simultâneas 


Brown Boveri 


Radiotelefones VHF/UHF 


para a transmissão de um canal telefónico e/ou sinais de 
telecomando 


Campos de aplicação: 

Correios, telégrafos, telefones 

Agências noticiárias 

Empresas de electricidade, de gás, de águas, 
petrolíferas 

Radiodifusão e televisão 

Policia, bombeiros e entidades oficiais 
Caminhos de ferro, aviação e navegação 
Industria, construção e transportes 


Representada em Portugal por: 
Sociedade de Electricidade 


À Groum Boveri Limitada 


Rua de Sá da Bandeira, 481-2º, Porto 


151660 a VI 


TÉCNICA XV 


TÉCNICA 


TÉCNICA XVI 


a técnica moderna emprega 


BETÕES REFRACTARIOS 


à base dos cimentos 


FONDU LAFARGE 
SECAR 250 


porque são 
práticos 
eficientes 
económicos 
REFRACTÁRIOS ATÉ 1800º C 
REFRACTÁRIOS ISOLANTES ATE 1700º C 
Os nossos serviços técnicos, especializados, estão à vossa inteira 


disposição para estudar, sem qualquer encargo da vossa parte, a 
adaptação dos BETÕES REFRACTÁRIOS a todos os vossos problemas 


AGUIAR & MELLO L.? 
P. do Município, 13-1.º — LISBOA — Tel. 321151/2 


NUMERO 423 


FEVEREIRO DE 1975 ANO L 


VOLUME XXXVI 


OBSERVAÇÕES SOBRE A 
«PROPOSTA PARA A REESTRUTURAÇÃO 
DOS CURSOS DE ENGENHARIA CIVIL» 


No N.º 420, de Setembro de 1974, desta revista, 
os professores do Curso de Engenharia Civil do |. 5. T., 
Eng.'(s) E. R. Arantes e Oliveira, A. Carvalho Quintela 
e À, F. Tovar de Lemos, publicaram uma proposta con- 
creta do elenco das cadeiras que constituirão o bacha- 
relato em Engenharia Civil. 


Foi com o maior interesse que lemos o programa 
apresentado, felicitando os autores pela sua precisão. 

Achamos contudo importante fazer alguns reparos 
ao modo como é proposto o ensino da cadeira de 
Materiais de Construção. Em primeiro lugar conside- 
ramos a Ciência dos Materiais de Construção como 
uma das cadeiras fundamentais da Engenharia Civil, 
a par da Resistência de Materiais. 


O engenheiro-civil não é um mero calculador de 
estruturas. No gabinete, ao projectar, tem de saber 
escolher e optar pelo material que melhor convenha 
para resistir a determinadas acções físicas e químicas. 
No estaleiro, ao construir, tem de conhecer os mate- 
riais de que se serve, 


Quando projecta uma estrutura, o engenheiro é obri- 
gado a ter noções precisas sobre o comportamento do 
material de construção à acção das forças. É portanto 
sob este aspecto que vai especialmente incidir o 
estudo dos materiais. Por esta razão, e por outras que 
referiremos mais à frente — a necessidade de fazer 
ensaios dos materiais — parece-nos indispensável que 
a cadeira seja dada depois de o aluno já possuir alguns 
conhecimentos de Resistência de Materiais. À sua situa- 
ção no 1.º ano não nos parece por isso adequada. 


Para a boa compreensão das propriedades fisicas, 
especialmente mecânicas, do material de construção, 
o aluno deve ter já estudado o ramo da Fisica que trata 
da constituição da matéria e da estrutura dos materiais, 
desde o nivel atómico ao macroscópico, a Física do 
estado sólido. 

Também pensamos que é desejável estudar as 
propriedades dos materiais de construção, propriedades 
essencialmente estatísticas, e em especial as do betão, 
depois de conhecer um pouco de cálculo estatístico, 
necessário para o controle da qualidade dos materiais. 
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A. DE SOUSA COUTINHO 
Investigador do Laboratório Nacional 


de Engenharia Civil 


O engenheiro civil precisa também possuir um 
certo conhecimento da tecnologia do fabrico dos mate- 
riais e muito especialmente da do betão, pois sendo o 
material de construção mais importante do nosso tempo, 
é obrigado a utilizá-lo, fabricando-o ou fiscalizando-o, 
praticamente, em todas as obras. 

Por isso deve o aluno ter noções suficientes de 
Quimica, que o habilitem.a conhecer os materiais metá- 
licos, e especialmente de química orgânica e de qui- 
mica dos silicatos, de modo que possa estudar conve- 
niegntemente os materiais plásticos, madeiras, tintas, 
betuminosos, etc. e os materiais siliciosos como a 
pecra, o solo, o cimento, o betão, os produtos cerã- 
micos, etc.. 

O conhecimento da Físico-química das superfícies 
é essencial para a compreensão dos fenômenos na super- 
ficia do grão dos materiais policristalinos, sejam ou 
não porosos ou microporosos, dos de rotura, de adsor- 
ção e desorção (tão importantes no comportamento do 
betão) e dos coloidais, que estão na base da estrutura de 
muitos materiais de construção. O seu estudo é também 
indispensável na Mecânica dos Solos. 

Com esta bagagem científica, que seria mais ou 
menos aprofundada conforme o grau do ensino — Bacha- 
relato, Licenciatura ou Doutoramente — o aluno poderia 
então estudar e compreender a segunda propriedade 
que o material de construção deve possuir em grau 
elevado — a sua durabilidade — ou seja a sua resis- 
tência à meteorização, englobando nesta a corrosão 
química e electroquimica, devida, quer à humidade, 
quer ao contacto com águas do mar, dos rios, dos 
solos, etc., as alternâncias de temperatura e humidade, 
as radiações solares ou de outra origem, as acções do 
metabolismo de certos agentes biológicos,etc.. 

Quanto à orientação a dar ao ensino da cadeira 
dos Materiais de Construção julgamos que teria de 
ser bem ponderada a sua estruturação, Até há umas 
dezenas de anos, digamos, até à primeira metade deste 
século, a Ciência dos Materiais de Construção consis- 
tia numa descrição das propriadades dos materiais; 
hoje tende-se para uma compreensão das propriedades 
em face da natureza da sua estrutura. 
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Estes dois aspectos não podem ser ignorados, 
Por isso torna-se indispensável que uma boa parte do 
ensino desta cadeira seja ministrado experimental- 
mente, no laboratório, Os alunos familiarizar-se-iam com 
as propriedades mais importantes dos materiais medin- 
do-as no laboratório; no caso do betão fabricá-io-iam 
mesmo para se aperceberem da sua natureza, tanto no 
estado fluido como no estado sólido, como aliás tem 
acontecido na Faculdade de Engenharia do Porto, 


Achamos portanto muito pouco o tempo que é 
destinado ao ensino dos Materiais de Construção, mesmo 
no bacharelato B meses — enquanto ao Desenho 
Técnico, que o engenheiro raramente tem de fazer na sua 
vida profissional, se dedicam três semestres. 


Compreendemos a dificuldade de introduzir no 
ensino da Fisica no nosso pais o ramo da estrutura da 
matéria e dos materiais, no sentido que é indispensá- 
vel não só ao engenheiro civil mas também ao mecânico 
e até mesmo ao electrotécnico. 

O mesmo se diz do estudo da Quimica e da Físico- 
-química das superfícies, cujos ramos mencionados 
julgamos não fazerem parte das matérias habitualmente 
ensinadas no nosso país nestas cadeiras de base, mas 
valeria a pena tentar introduzi-los. 

Por isso, e enquanto estas disciplinas não forem 
leccionadas como foi preconizado — o que implica, 


como dizem os autores do artigo em discussão «o reco- 
nhecimento de que a maioria dos engenheiros exerce 
a sua profissão sem recurso à formação teórica» (em 
vez de teórica deve ler-se cientifica básica) «que a escola 
pretendeu dar-lhes» — é indispensável que o professor 
de Materiais de Construção transmita aos alunos tais 
conhecimentos basilares, de uma forma que terá de 
ser muito condensada. Isto agravará mais ainda o 
tempo que é necessário dispensar a esta cadeira. 


Sabemos bem das dificuldades do ensino de 
Materiais de Construção tal como foi apontado e com 
larga participação de aulas práticas, de laboratório. 
A inexistência de um ensino experimental no nosso 
País torna efectivamente difícil, mas não impraticável, 
um tal esquema. 

Em resumo, a disciplina da Ciência dos Materiais 
de Construção deve ser ministrada depois de o aluno 
possuir já conhecimentos adequados de Fisica, Quimica 
e especialmente de Resistência de Materiais e abarcar 
pelo menos três semestres, pois não faz sentido que 
no elenco proposto a Resistência de Materiais e a 
Análise das Estruturas abranjam 4 semestres, o Desenho 
3 semestres e os Materiais de Construção, disciplina 
que, voltamos a insitir, consideramos também funda- 
mental na Engenharia Civil, apenas um semestre, 
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Feira Internacional de Plásticos + Caucho 


Feira Técnica Industrial 


Salão das Industrias de Equipamento e Manutenção Automóvel Paris 


World Congress ou Finite Element Methods in Structural Me- 
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43.º Congresso Internacional de Fundição 
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| LOCAL, DATA, ORGANIZAÇÕES 
E INFORMAÇÕES 

| Dusseldorf 7 - 14 Outubre --75 

| Rep. Fed. d'Allemagne 

| 


Tel - Aviv 
Outubro 


2a 19 de Outubro - 75 


Bomnemouth, Dorset, England 121th - 171th 
October - 75 at the Royal Bath Hotel 


Dússeldort 22 al 29 Mayo - 76 
Minéria 76 

Basle - (Suiça) 
22 a 26 Junho - 76 


Bucareste 
5 a 10 Setembro - 76 
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Resumo dos artigos publicados na «Técnica» nº 423 


Ano L — Fevereiro 1975 


C. D.U. 622.003 
C. DINIS DA GAMA 


BASES DO ESTUDO ECONÔMICO 
DE UM EMPREENDIMENTO MINEIRO 


Técnica N.º 423 — L— 2-1975, p. 111 a 120 


À análise econômica aplicada a empreendimentos mineiros 
pressupõe o domínio de uma sequência de conhecimentos 
interligados, respeitante a informação, avaliação e decisão, 
à validade daquela análise é lógicamente condicionada 
pela representabilidade destas actividades. Neste trabalho, 
faz-se referência à necessidade de obtenção de dados de 
confiança, mensionam-se os vários métodos determinísticos 
e probabilisticos de avaliação econômica e apontam-se 
critérios de dicisão mais utilizados, com especial rea'ce 
para os que se apoiam na simulação por computador. 

A natureza interdisciplinar deste estudo é posta em relevo, 
embora certos aspectos da análise sejam típicos da Indústria 
Mineira. 


C. D. U, 624.022,327.002.22:699,841 


5. FERRY BORGES, A. SARAIVA, 
J. PEREIRA E V. MONTEIRO 


METODOLOGIA PARA O ESTUDO SISMICO 
DE EDIFICIOS PREFABRICADOS POR PAINEIS 


Técnica N.º 423— L— 214975, p. 121 a 126 


à homologação de edifícios prefabricados por painéis em 
zonas sismicas requer a realização de vários tipos de 
estudos conducentes à caracterização do comportamento 
das estruturas dos edifícios sob acções sismicas, Tais 
estudos dizem respeito principalmente a: i) idealização 
do comportamento estrutural sob a acção de forças hori- 
zontais; ii) análise dinâmica em regime linear e não-linear; 
Hi) resisitência de ductilidade dos elementos estruturais 
isolados ou da estrutura global. No tratamento desses 
problemas torna-se indispensável conjugar de forma ade- 
quada a utilização de técnicas tanto analíticas como 
experimentais, 

Nos últimos dez anos têm sido realizados pelo Labo- 
ratório Nacional de Engenharia Civil (LNEC) diversos 
estudos sobre o comportamento sísmico de edifícios pre- 
fabricados por painéis, a construir em Portugal e no estran- 
geiro, Sintetizam-se no presente trabalho os principais 
aspectos da experiência adquirida ieste domínio e a meta- 
dologia que tem sido seguida. 


C. D. U. 518.6 
ALVARO BENTO LEAL 


SOLUÇÃO DA EQUAÇÃO DE POISSON PELO METODO DE 
RITZ-PROCESSO GERAL PARA QUALQUER CONTORNO 


Técnica N.º 423 — L-— 24975, p. 127 a 132 


É aplicado à equação de Poisson bidimensional com condi- 
ções na fronteira de Dirichlet a técnica da semi-discreti- 
zação. Mostra que o sistema de equações diferenciais 
obtido permite uma formulação variacional equivalente. 
à solução do problema variacional pelo método de Ritz, 
reduz-se ao cálculo de um conjunto de funções que se 
exprimem explicitamente em termos de coordenadas de 
pontos sobre o contorno que limita o domínio da solução. 


Cc. D. U. 621.317.34 
J. M. VIEGAS GONCALVES 


IDENTIFICAÇÃO DE SISTEMAS LINEARES 
POR AMOSTRAGEM FINITA NO DOMINIO DO TEMPO 


Tócnica N.º 423— L— 2.1975, p. 133 a 140 


Apresenta-se um método para identificação de funções 
caracteristicas de sistemas descritos por equações diferen- 
ciais ordinárias lineares de coeficientes constantes a partir 
dum número finito de valores da resposta impulsiva toma- 
dos em pontos equidistantes no domínio do tempo. 
Indica-se a estrutura dum programa para o cálculo auto- 
mático dos polos e zeros da função característica a 
identificar a partir dos valores observados. 

O método é aplicado à determinação experimental da 
função de transferência dum filtro LCR passa-baixo de 
ordem 5. 


C. D. U. 551.501.::629.19:681.14 
RODRIGO MARTINS E ISAURA MARTINS 
CONVERSOR DIGITAL ANALOGO TEMPORAL 
Técnica Nº 423 — L— 2.1975, p. 141 a 164 
Ensala-se a síntese de um pequeno computador especial- 
mente projectado para controlar os movimentos de uma 


antena helicoidal, destinada à captação de fotografias 
meteorológicas enviadas pelo satélite n«ESSA - 8». 
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U. D. C. 621.317.34 
J. M. VIEGAS GONÇALVES 


LINEAR SYSTEM IDENTIFICATION BY FINIT SAMPLING 
IN TIME DOMAIN 


Técnica N.º 423— L— 2.1975, p. 133 a 140 


This paper presents a method for characteristic function 
identification in systems which are described by sets 
of linear ordinary differential equations with constant 
coefficients, based on finit sampling in time domain of 
tha system's impulse response. 

The setructure of a program used to compute the location 
of poles and zeros from the sampled values is described. 
The method is applied to the experimental determination 
of a fifth order LCR low pass filter transfer function. 


U. D. €. 551.501:629.19:681.14 
RODRIGO MARTINS e ISAURA MARTINS 
ANALOG TO TIME DIGITAL CONVERTER 
Técnica No 423 — L— 2.1975, p. 141 a 164 


«lt is diseussed the design of a special purpose computer, 
capabla to control two movements of an helycoidal antena 


which receives meteorological photographies from ESSA - 8 


satelyte»., 


U. D. €. 622.003 
C. DINIS DA GAMA 


BASIS FOR THE ECONOMIC ANALYSIS 
OF A MINING VENTURE 


Técnica No 423— L— 24975, p. 111 a 120 


Tha economic analysis of mining prospects is conducted 
through the knowledge of a series of topics concerning 
information,evaluation and decision. The validity of each 
analysis depends upon the reliability of those forecasting 
activities. 

In this paper, the need for obtaining reliable data, as well 
as for apolying the available methods, both deterministic 
and probalilistic, of economic appraisal, and also to utilize 
decision criteria are underlined. Special emphasis is 
ascribed to computer simulation methods. 

The interdisciplinary nature of this economic analysis is 
mentioned, although some aspecis are characteristic of the 
Mining Industry. 
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METHODOLOGY FOR SEISMIC STUDIES 
OF PREFABRICATED PANEL BUILDINGS 


Técnica Nº 423— L— 24975, p. 121 a 126 


Tho agreement for the use of prefabricated panel buildings 
in seismic areas requires different kinds of studies in order 
to determine the major aspects of the seismic behaviour 
of the buildings. Such studies should concern: i) idealization 
of structural behaviour under horizontal loads; Ii) linear 
nud non-linear seismic response; and iii) resistance and 
ductility of individual members and the whole structure. 
Suitable combination of analytical and experimental techni- 
ques should be used for studying these problems. 

In the last ten years the Laboratório Nacional de Enge- 
nharia Cívil (LNEC) has studied several prefabricated 
panel structures to be built both in Portugal and abroad, This 
paper briefly reports the experience of the LNEC in this 
fiald and describes the methodology that is being used. 
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SOLUTION OF POISSON EQUATION BY RITZ METHOD. 
GENERAL PROCEDURE FOR ANY BOUNDARY 


Técnica N.º 423— L— 21975, p. 127 a 132 


The straight line method is applied to bidimensional Poisson 
equition with Dirichet boundary conditions. |t shows that 
the resulting system of differential equations allwos an 
equivalent variational statement. The solution of the varia- 
tional problem by Ritz method is reduced to the calculation 
of a set of function which are closed forms in terms of 
tha boundary point coordinates. 
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